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Abstract 

We derive general bounds on operator dimensions, central charges, and OPE coeffi- 
cients in 4D conformal and N = 1 superconformal field theories. In any CFT containing 
a scalar primary (ft of dimension d we show that crossing symmetry of {4>4>4>4>) implies 
a completely general lower bound on the central charge c > f c (d). Similarly, in CFTs 
containing a complex scalar charged under global symmetries, we bound a combination 
of symmetry current two-point function coefficients t ij and flavor charges. We extend 
these bounds to J\f = 1 superconformal theories by deriving the superconformal block 
expansions for four-point functions of a chiral superfield $ and its conjugate. In this 
case we derive bounds on the OPE coefficients of scalar operators appearing in the 
$ x & OPE, and show that there is an upper bound on the dimension of when 
dim<I> is close to 1. We also present even more stringent bounds on c and r . In 
supersymmetric gauge theories believed to flow to superconformal fixed points one can 
use anomaly matching to explicitly check whether these bounds are satisfied. 
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1 Introduction 



Near-conformal dynamics may describe physics beyond the Standard Model, in addition 
to facets of QCD itself. Examples include walking [1-6] and conformal [7-9] technicolor, 
dynamical explanations of the flavor hierarchies [10-13], solutions to the SUSY flavor prob- 
lem [14-27], solutions to the \ijB\i problem [28-32], and so on. While many of these ideas 
are promising, they often rely crucially on assumptions about the behavior of strongly- 
coupled field theories. However, conformal symmetry itself severely restricts the structure 
of these theories, and it is not fully understood which assumptions are consistent with these 
restrictions and which are not. 

In [33-35] significant progress was made in understanding the range of behavior that is 
possible in 4D conformal field theories. The key insight is that crossing symmetry of four- 
point functions requires that coefficients appearing in the operator product expansion (OPE) 
not be too large. Combined with certain assumptions about the spectrum of operators, these 
constraints can potentially lead to a contradiction with unitarity, allowing one to rule out 
the spectrum. Concretely, in [33, 34] it was shown that there is a completely general upper 
bound on the dimension of the lowest- dimension scalar primary operator appearing in the 
OPE x of a real scalar primary of dimension d with itself, < f(d), where f(d) is a 
function that is determined numerically. In [35] it was also shown that one could compute 
an upper bound on the coefficient of the three-point function (<fxpO) for any scalar primary 
O appearing in the OPE. 

In the present work, we extend the analysis of [33-35] in several directions. First, we 
examine crossing symmetries of correlators involving charged fields in CFTs with global U(l) 
symmetries, focusing in particular on the additional constraints that are present in supercon- 
formal theories. We consider a chiral superconformal primary operator <3> of dimension d, and 
show how the four-point function may be expanded in terms of "superconformal 

blocks", which sum up the contributions of a given superconformal multiplet appearing in 
the $ x $' OPE. Since each superconformal multiplet contains a finite number of primary 
operators under the conformal sub-algebra, superconformal blocks may be decomposed into 
a finite sum of conformal blocks. While such a decomposition was previously known in the 
context of M = 2 and M = 4 theories [37], we believe that the M = 1 result we present 
is new. We further show how the M = 2 superconformal blocks derived in [37] may be 
decomposed in terms of M = 1 superconformal blocks, providing a non-trivial check on our 
result. 

Second, we combine our superconformal block analysis with the methods of [33, 34] to 
derive bounds on the spectrum of operators appearing in the <3> x <3>t OPE. In particular, 
we find that there is an upper bound on the dimension of the $t<i> operator (defined as 
the lowest-dimension scalar appearing in $ x $t) when d is close to 1. Since the chiral 
operator $ 2 with dimension 2d always appears in $ x $, one cannot reproduce our bound 
on <£>t$ by simply applying results from [33, 34] to the real or imaginary parts of $. We also 
compute bounds on the OPE coefficient of any scalar superconformal primary appearing in 
$ x $t, independent of assumptions about the spectrum. Our dimension and OPE bounds 
constitute completely general non-perturbative results about non-BPS quantities in Af = 1 
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superconformal theories. 

Third, we use crossing relations among complex scalars to study OPEs involving con- 
served currents, both in the supersymmetric and non-supersymmetric context. If is a 
complex scalar primary of dimension d, then the OPE x 0* contains global symmetry 
currents J aI , whose coefficients are fixed by a Ward identity to be proportional to the 
charges of 0. Crossing symmetry of (00*00*) then implies an upper bound on the charges of 
relative to the "flavor central charges" t ij , defined by the coefficient of the two-point function 
(J aI J bJ ) oc t ij . Specifically, we show tjjT i T j < f T (d), where tjj is the inverse of r IJ , and 
T 1 are the global symmetry generators in the representation. We further strengthen this 
bound when is the lowest component of a chiral multiplet $ in a superconformal theory, 
in which case flavor currents appear as descendants of scalar operators J 1 . 

Last, we turn to OPE's involving the stress tensor T ab in both supersymmetric and non- 
supersymmetric theories. We show that in any CFT containing a real scalar primary operator 
of dimension d, there is a completely general lower bound on the value of the central charge 
c > f c (d). This again occurs because crossing symmetry of the four-point function (0000) 
requires that the OPE coefficient in front of the stress tensor 0(x)0(O) ~ T ab (0) not be 
too large, and this coefficient is fixed in terms of c and d by a Ward identity. This can 
perhaps be viewed as a four-dimensional counterpart to the bound on c derived in [36] for 
two-dimensional CFTs. Once again, we strengthen this bound in the supersymmetric case, 
where the stress tensor appears as a descendant of the U(1)r current in the $ x $^ OPE. 

Our bounds on tijT i T j and c are particularly interesting in supersymmetric theories 
since these quantities are determined in terms of the superconformal U(1)r symmetry as 
r IJ = -3Tr( J RT 7 T J ) and c = i (9Tri? 3 - 5Tri?), and may be calculated via 't Hooft 
anomaly matching. Since <3> is chiral, its dimension is also determined in terms of the U(1)r 
symmetry as d = Thus, one may check whether these bounds are satisfied in the myriad 
asymptotically free M = 1 theories that are believed to flow to superconformal fixed points, 
and we will demonstrate that this is the case in a few simple examples. 



2 Preliminaries 

2.1 CFT Review 

We will begin our discussion by reviewing some basic facts about 4D conformal field theories. 
The conformal algebra may be written as 

[M ab , P c ] = P a r] bc - P b r] ac , [M ab , K c ] = K a r] bc - K b r] ac 

[M ab , M cd ] = r] bc M ad - r] ac M bd - r] bd M ac + rj ad M bc 

[D,P a ] = P a , [D, K a ] = -K a 

[K aj P b ] = 2 Vab D-2M ab , (2.1) 
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and primary operators (9 7 (0) are defined by the condition i^" a (9 7 (0) = 0, 1 where K a is 
the generator of special conformal transformations. Fields may then be constructed by 
exponentiating the translation operator, 0\x) = e xP O I (0). Here I denotes possible Lorentz 
indices, which can be labeled by according to the representation of SO(4) = SU(2) x 
SU(2). For example, traceless symmetric tensors O ai - a '(x) have j = j = 1/2 (which we 
call the "spin-/" representation). We refer the reader to Appendix A for a more complete 
summary of the conventions used in this paper. 

In 4D conformal field theories the correlation functions of primary operators are highly 
constrained (see e.g. [38]). In particular, the two-point function for a spin-/ primary operator 
O ai " Ml (x) of dimension A can in general be written as 

I aibl (x 12 ) . ..I a ' bl (x 12 ) 

x 12 

rf - 2^-, (2.2) 

x z 

where xu = x\— x 2 , and the indices a± . . . ai and b\ . . . b\ are implicitly symmetrized and made 
traceless. Unitarity requires that the coefficient of the two-point function is positive, so that 
one can choose a basis of primary operators with the above normalization, where additionally 
two-point functions between different basis elements are taken to vanish, (0(x-y)0 '(£2)) — 
for O 7^ O' . Positivity of the two-point functions of descendant operators then further 
imposes the unitarity bounds [39] 

A > 1 (1 = 0), 

A > 1 + 2 (/>!)• (2.3) 



(O ai - a '(x 1 )O bl - bt (x 2 )) = 
I ab (x) = 



Three-point functions between scalar primary operators 4>i(x) of equal dimension d and 
a spin-/ primary O ai " Ml (x) of dimension A are fixed up to an overall constant as 

z a 

If we take <fti = <p 2 , it is straightforward to see that invariance under x\ -H- x 2 requires that / 
must be even in order for the three-point function to be non-vanishing. However, if cf>i 7^ <p 2 
then odd-/ primaries are also allowed. Note that Lorentz representations with j 7^ j cannot 
appear because there does not exist a function built out of the Xj's having the required 
transformation properties. As reviewed in [33], when <pi and <p 2 are real, the coefficients 
^4>i<ho are necessarily real in the basis of Eq. (2.2). 

1 For notational convenience we will leave the adjoint action of K a implicit in expressions such as this, so 
that -ftT a O / (0) — > [K a , O / (0)]. Fcrmionic gradings should be respected. For example if O is bosonic, then 
Q 2 is short for { Q & , \Q & , O] } . 



A</, 1(j 
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Finally, four-point functions of scalar operators are not completely determined by sym- 
metry considerations alone, and in the case of equal dimensions can always be written as 



glu, v) 

(0l(^l)02(^2)^3(^3)04(^4)) = od^d ' 

(2.5) 

X 12 X 34 

2 2 2 2 

where g(u, v) is a function of the conformally-invariant cross ratios u = x l 2X % 4 and v = x ^ x f . 

Though g(u,v) is not fixed by conformal symmetry, it is fully determined by the dynamical 
data of the theory, namely the spectrum of operator dimensions and spins A, I and three- 
point function coefficients X^^o- This is most easily seen through the operator product 
expansion (OPE), which relates a product of operators at different positions to a sum over 
operators at a single position. In the case of scalar primaries, we can write 

<Mx)0 2 (o) = Yl C/(z,P)e>'(o), (2.6) 

where I stands for possible Lorentz indices. We use the notation O £ <p\ x <f) 2 to mean that 
the sum should be taken over primary operators occurring in the OPE of 0i with <p2- The 
operator Ci(x, P) may for example be determined by inserting the OPE into the three-point 
functions and using the known form of the two-point functions [40]. 

Taking the 0i(xi) x ^2(^2) OPE and the 03(0:3) x ^4(^4) OPE in the four-point function 
then leads to the conformal block decomposition 

g(u,v) = ^2 X^ 2 oX^409a,i( u ^ v )^ ( 2 - 7 ) 

0£<j>lX<t>2 

where the "conformal blocks" gA,i{u,v) are given explicitly by [40] 

(— 1)' zz 

9a,i( u > v ) = — 77i = [ k A+i(z)k A _i_ 2 (z) - z z] 

A 1 z — z 

kp{x) = xP/ 2 2 F 1 ((3/2, (3/2, /3; x), (2.8) 

and the change of variables u = z~z and v = (1 — z)(l — z) has been used. We note in passing 
that the conformal blocks can also be elegantly derived by viewing them as eigenfunctions 
of the quadratic casimir of the conformal group [41]. 

If we take all of the scalars to be identical, then invariance of Eq. (2.5) under X\ ■£> 23 
leads to the "crossing symmetry" constraint 

Y X l<t>o9A,i( u , v ) = (^) Y X 1<I>o9a,i(v,u), ( 2 - 9 ) 

which must satisfied by any consistent spectrum of dimensions, spins, and choice of three- 
point function coefficients. A key point is that unitarity requires X^q £ so the coefficients 
X 4><t>o appearing above are positive. Invariance under xi -H- x 2 again tells us that only even- 
spin operators may appear, and other exchanges do not give any new information. 
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2.2 Bounds from Crossing Relations 



In [33, 34], the crossing relation of Eq. (2.9) was used to derive an upper bound on the 
dimension of the lowest-dimension scalar operator appearing in the OPE 0x0. In [35] 
bounds were also derived on the size of the three-point function coefficients of scalar operators 
appearing in x 0. The techniques employed depend on the explicit expression Eq. (2.8) for 
conformal blocks, together with the unitarity requirement A;L, > 0. We now review these 
techniques; in the following subsection we will discuss some generalizations. 

Let us begin by showing how to bound the OPE coefficient-squared \ 2 0q = A|^ 0o of a 
given operator Oq of dimension A and spin l appearing in x 0. We first rewrite the 
crossing relation by separating out and dividing by the contribution of the unit operator, as 
well as separating out the contribution of the particular operator Oq whose OPE coefficient 
we would like to study, 

X 2 o F Ao , lo (u,v) = 1- X 2 oF A>l (u,v), (2.10) 

where 

F A ,i{u,v) = , (2.11) 

U — V 

and we have used that g o(u,v) = 1 for the unit operator. Note that F A ^ depends on d, 
though we are suppressing this dependence for brevity. Eq. (2.10) is a linear equation in the 
space V of functions of two variables which are invariant under u ■<->• v. It encodes an infinite 
number of relations between OPE coefficients A©, but general statements about solutions 
consistent with unitarity (A|, > 0) can be difficult to extract. The approach of [33-35] is to 
consider a real linear functional a G V* = Hom(V,M), which satisfies 

a(*A ,Jo) = 1, and ( 2 - 12 ) 
ot(F A j) > 0, for all other operators in the spectrum. (2-13) 

Then applying a to both sides of Eq. (2.10), we obtain a bound 

\ 2 O0 = a(l) - A 2 a(F Aji ) < a(l), (2.14) 

where we have used that \ 2 > by unitarity. Let us denote by S the subspace of a G V* 
which satisfy the constraints (2.12, 2.13) (depicted in Figure 1). In many cases of interest, 
S is non-empty, so a non-trivial bound on the OPE coefficient-squared \@ exists. 

Given bounds on Xq o , bounds on the dimension of Oo may or may not follow as a 
consequence. For example, suppose we assume that Co with dimension A is the lowest- 
dimension scalar appearing in x 0. Then if we can find some a such that \ 2 0o < a(l) < 0, 
then we have found a contradiction with unitarity, implying that it's impossible that Oq has 
dimension A . 
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Now our bound \@ < is most interesting when a(l) is as small as possible. Thus we 
would like to minimize a(l) over all a G S. This problem resembles an infinite-dimensional 
version of a linear program, which usually refers to a linear optimization problem over M. n , 
subject to a finite number of affine constraints. Linear programs have been well-studied in 
mathematics and computer science, and a number of efficient algorithms for their solution are 
known. A key observation is that since the search space is an intersection of half-spaces (one 
for each inequality) and hyperplanes (one for each equality), it is convex. Consequently, the 
optimum of any linear function lies on the boundary of the search space, and can be reached 
deterministically by following the direction of steepest descent (either along the boundary 
or in the interior). 




Figure 1: The "search space" S C V* is the intersection of the hyperplane oi(Fa ,i ) = 1 with the 
convex cone of linear functionals a satisfying a(Fj\ i) > for all (A, I) in the spectrum. 

A first step towards making our problem tractable via these methods is to restrict to 
a finite-dimensional subspace W C V*. Then, minimizing a(l) over a G W D S will give 
a possibly sub-optimal, but still valid bound \ 2 Qq < a(l). The choice of W is somewhat 
arbitrary and unfortunately can have a significant effect on the answer. A convenient class 
of subspaces is given by taking linear combinations of derivatives at some point in z, z space. 
Following [33-35], we take these derivatives around the point z — z — 1/2 (which is invariant 
under u •f-)- v). That is, we define W4 C V* to be the space of functionals 

a:F(z,z) m- J2 a mn d^F(l/2, 1/2) (2.15) 

m+n<2k 

with real coefficients a mn . 2 We can then scan over M4 PI S by varying the a mn , subject to 
the constraints of Eqs. (2.12, 2.13). One hopes that as we take k — > oo, our search will cover 
more and more of S, and our bound will converge to the optimal one. 3 

2 In addition to being simple to describe, the spaces Wk are computationally convenient, since one can in 
fact derive relatively simple analytic expressions for derivatives of the functions Fa,i{z,z) at z = z = 1/2, 
and these expressions can be computed efficiently using recursion relations (see Appendix B). 

3 Hcre, "optimal" means "optimal given our assumptions," namely the diagonal crossing relation Eq. (2.9) 
and unitarity of each OPE coefficient in 0x0. These are a small subset of the full consistency relations of 
a CFT, so it's certainly possible that inputting more information could lead to even stronger bounds. 
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Even after restricting to Wfe, our problem differs from a typical linear program in that 
Eq. (2.13) includes an infinite number of affine constraints on a. For example, if we are 
interested in bounding \ 2 0o with no additional assumptions on the spectrum, then we must 
demand a(F^i) > for all (A, I) obeying the unitarity bound. Alternatively, if we wish to 
bound the OPE coefficient of the lowest-dimension scalar in <fi x <ft, we must take a(F^i) > 
for all scalars with A > A , and all (A, I) with I > that obey unitarity. In each case, we 
have a continuously infinite number of constraints on a — one for each (A, /) pair. 

In a typical linear program, the search space is a convex polytope in R n , given by an 
intersection of a finite number of half-spaces and hyperplanes. In our case, the search space 
S C V* is still convex, since it is an intersection of half-spaces Ua,i = {« : Q;(Fa,j) > 0} 
and a hyperplane H = {a : <y(F& 0i i ) = 1}. However, because A can vary continuously, 
S is not a polytope. In general, the intersection S D W with any finite-dimensional W is 
"piecewise-curved," e.g. it has (not-necessarily flat) faces, (possibly curved) edges, vertices, 
etc. Consequently, we expect that at finite k, as we vary the underlying parameters of our 
problem (d, A , etc.), our bound will vary in a "piecewise-curved" way, with corners as the 
optimal a = a* passes over edges on the boundary of S n Wt- 

In order to apply linear programming techniques, we need to approximate S fl V\4 by a 
polytope. 4 Let us pick some finite discrete set D = {(Aj,/j)} and reduce the constraints in 
Eq. (2.13) to simply a^F^^) > for all (Aj, U) G D. This expands the search space, and we 
are now in danger of obtaining an invalid bound if the optimal a = a* satisfies a*(i<A'/) < 
for some (A', /') not in D. However, this danger disappears as we increase the size of D 
and approximate S fl W4 by more and more refined polytopes. A type of discretization that 
works well in practice is 

D = {(A min + ne,Z) : n = 0, . . . , iV and Z = 0, 2, . . . , L}, (2.16) 

where iVe and L are large numbers (say ~ 50), and e is some small step size (say e ~ .1 or .01). 
By decreasing e, we can ensure that violations of our constraints a*(i 7 A',« / ) < become less 
and less important. One must also ensure that a{F^ {) is greater than zero asymptotically 
as A, / — > oo. This is easy to check using the analytic expressions for derivatives of F^j 
given in Appendix B. In practice, the optimal a = a* often obeys the asymptotic constraint 
automatically, provided iVe and L are sufficiently large. 

2.2.1 Solutions to the Crossing Relations from Linear Programs 

In this subsection, we will show how in principle, the linear program described above produces 
not just an OPE coefficient bound, but also the corresponding "optimal" solution to the 
crossing relations consistent with the given assumptions. This type of solution doesn't 
necessarily have anything to do with CFTs, since we're only inputting a subset of the full 
CFT consistency relations. However, we mention it here because it helps give some intuition 

4 Actually, there do exist algorithms to solve more general classes of "convex optimization problems", 
which can involve curved search spaces. It might be interesting to investigate whether any of these can be 
applied to crossing relations. 
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for properties of optimal solutions a* G S. The results of this subsection are not used 
elsewhere in the paper, so the reader should feel free to skip to Section 2.3 if desired. 

Let us briefly introduce some notation. A subset K C V of a finite-dimensional real vector 
space V is called a convex cone if A1X1+A2X2 G K for all Xi, X2 G K and Ai, A2 G K+. The dual 
cone of is the space of linear functions K v = {£ G V* such that £(%) > for all x G if}. 
One can show that if K is closed, then (if v ) v = K. 

In the following, let us be cavalier and pretend that the space V* of possible en's is finite- 
dimensional. Suppose we have run our linear program and arrived at an optimal a* G S 
that minimizes a(l) subject to Eqs. (2.12) and (2.13). Since a* lies on the boundary of 
S, some of the constraints defining S must be saturated at a*. That is, there is some set 
{FAi,k for % — 1, 2, . . . } such that a*{F Au i^) = 0. We also have a*(F Aot i ) = 1 by assumption. 

For the local geometry of S near a*, only the constraints which are saturated at a* are 
important. In other words, we can imagine locally replacing S with the set a* + /C a „, where 
/C Q< , is the convex cone 

£ a , = {/3 : f3(F AiM ) > for % > 1, and /3(F Aoilo ) = 0}. (2.17) 

Note that /Cq,, is the dual cone of 

X£ = RF Ao , io + ^R + F AsA , (2.18) 

namely the positive span of the F/^.^, plus -Fa ,Zo with an arbitrary real coefficient. 

Now, the condition that a* minimize a*(l) means that (a* + 5a) (1) > a*(l) for all 5a 
pointing into the interior of S, that is all 5a G IC at . But this just means that 1 is an element 
of the dual cone K,^ , so that there exist coefficients q G M and Pi G IR+ with 

1 = qF Ao ,i + J2Pi F ^U- (2-19) 

In other words, the (Aj,/j) whose constraints are saturated at a*, along with (A ,/o) ; gi ye 
the spectrum of a solution to the crossing relation Eq. (2.10). Further, if q > then this 
solution is consistent with unitarity, since the Pi are positive. 

A useful geometric picture (illustrated in Figure 2) for arriving at the above result is to 
imagine picking a metric and thinking of —F Aij i. as specifying "normal vectors" to the search 
space S at a*. The minimum of a(l) occurs precisely when the vector — 1 G V (which points 
in the direction we want to go) is in the positive span of normal vectors to S. Meanwhile, 
since we can never move off the hyperplane a(F A(h i ) = 1, it doesn't matter whether —1 has 
a component in the direction of F Ao j , which is why q can have either sign. 

In practice, we must solve our linear program by first restricting to a finite dimensional 
search space S R Wfe. In this case, the optimal a* will have a few saturated constraints 
a*(F Au [.) = (i — 1, . . . , N k ). We can see this explicitly in Figure 3 which plots a*(F A j) for 
varying A and /, where a„ is the solution to a linear program. Note that a*(F Aj[ ) has zeros 
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Figure 2: Picking a metric on V*, we arrive at the following picture. The linear functional that 
minimizes a(l) is the unique point a* on the boundary of S where —1 is in the positive span of 
the "normal vectors" —F^.j. to S at a*. Here, the three parallel arrows illustrate the direction 
of steepest descent of a(l). We have suppressed an infinite number of dimensions (including the 
F& j direction) in order to draw this figure in the plane. 

at particular (A, I), but of course never becomes negative. As we increase k, we expect new 
zeros of a* (-Fa,;) to appear, with Nk eventually running off to infinity. If the limit k — > oo is 
"wcll-bchavcd" in some appropriate sense, we might hope that the zeros (A,. /,) for small A 
and I values converge quickly as k — > oo, giving us some information about the low-dimension 
and spin part of the spectrum corresponding to the "optimal" solution Eq. (2.19). Indeed, 
this seems to be the case in practice. It would be interesting to see if this information has 
any practical applications. 

2.3 Limitations and Generalizations 

One limitation of the formalism outlined above and the one used in [33-35] is that one 
only learns about the OPE of a real scalar with itself. In particular, the formalism does 
not allow one to distinguish between operators appearing in the OPE that have different 
global symmetry charges. For example, in TV" = 1 superconformal field theories there is a 
global U{1)r symmetry, and chiral operators have dimension d = If we take to be 
the lowest component of a chiral multiplet $, then the Re[0] x Re[0] OPE will contain both 
operators in the x OPE having U(1)r charge and operators in the x 0* OPE 

that are neutral under U(1)r. Since the 2 operator appearing in the 0x0 OPE is chiral, 
it always has dimension 2d and automatically satisfies the bounds derived in [33, 34]. Thus, 
we unfortunately do not learn anything new about the £/(l)R-singlet non-chiral operators 
appearing in x 0*. 5 

With this motivation in mind, let us consider more carefully what crossing relations 

5 Another example, extensively discussed in [33], is that the bounds do not distinguish between SU(2)- 
singlct and S'[/(2)-triplet operators appearing in the h x h OPE in conformal technicolor scenarios. 
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Figure 3: A plot of a*(F A i) for various A,/, where a* S S gives the strongest bound on the 
OPE coefficient of the lowest-dimension scalar Oq G <j) x <j>. Here, we have taken dim</> = 1.1, 
Ao = dim Co — 115, and k = 4. We show only pairs A, I satisfying unitarity. Note that at*(F& i) 
is never negative in this range, consistent with the constraints of our linear program, although it 
has zeros (A,/) G {(2.2,0), (8.2,0), (4,2), (7.7,2), (6,4)}. Note also that a*(F Ao , ) = 1, as required. 



apply in the case of a complex scalar charged under a global U(l) symmetry. We must 
first determine the unitarity constraints for three-point functions involving a complex scalar. 
We know that the correlator with a real spin-/ operator (<j)(xi)4>*(x2)O ai - ai (xs)) must be 
invariant under the exchange X\ -H- £2 combined with complex conjugation. From this we 
learn that even-/ operators must have coefficients that are real, \^*o = A^„ , and odd-/ 
operators must have coefficients that are imaginary, X^o = On the other hand, 

the three-point function ((f)(xi)(f)(x2)O ai '" ai *(x 3 )) must simply be invariant under x\ <H- 22, 
and hence only even-/ operators may appear. In this case, however, the coefficient \^o* 
is in general complex. Of course, the above arguments reproduce exactly what we would 
have concluded by breaking <fi into its real and imaginary parts <fi = <fii + ifo, and using the 
requirement from Section 2.1 that X^ o £ 

Now let us consider the four-point function {(f)(xi)(f)(x2)4>* (x 3 )4>* (x 4 )) . We can evaluate 
this in two qualitatively different ways: by taking the (j>{x\) x 0(x 2 ) and 0*(x 3 ) x 0*(x 4 ) 
OPEs, or alternatively by taking 0(xi) x 0*(x 4 ) and 0(x 2 ) x 0*(x 3 ). Equating the resulting 
expressions leads to the crossing relation 

^2 \^UO*\ 2 9a,i{ u i v ) = (-) \ X Wo\ 2 9a,i( v ^ u )- ( 2 - 20 ) 

Although one could conceivably apply the ideas of Section 2.2 to this kind of relation, we 
have had more success applying linear programs to crossing relations that display symmetry 
under u -H- v — that is, relations which involve the same spectrum of operators on both sides. 
As we will see, Eq. (2.20) implies two such independent relations that must be satisfied in a 
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consistent theory. By adding the equation to itself we can immediately derive one of them, 

\ x o\ 2 9A,i(u>,v) = (^) Yl M 2 gA,i{v,u), (2.21) 

0£<f>X(f> 0£(f>X(j> 
OG?ix0* 0£(px<f>* 

where Xq is shorthand for the appropriate three-point function coefficient. 

The simplest way to see the second crossing symmetry constraint is to alternatively relabel 
the coordinates and consider expanding the four-point function {(j)(xi)4>*(x2)<j)(x 3 )(j)*(x4 : )) by 
taking the 4>{x\) x 0*(x2) OPE and the 0(0:3) x <f)*(x4) OPE. Exchanging x\ <-> £3 then leads 
to the constraint 

\^*o\ 2 {-l) l gA,i{u,v) = \^4>*o\ 2 (-l) l gA,i{v,u), (2.22) 

which is a crossing symmetry relation that only involves operators in the 0x0* OPE. 
Here the (—1)' factors appear because the odd-/ operators necessarily have 3-point function 
coefficients that are imaginary, and their square (and not absolute value squared) enters the 
conformal block decomposition when the 0's have the above ordering. One can show that 
Eq. (2.22) in fact follows from Eq. (2.20) through repeated use of the identity gA,i{ u , v ) — 
(~^) l gA,i{u/v, 1/v) along with the knowledge that only even-/ operators appear in the x 
OPE. It can similarly be verified that other crossings do not contain any additional 
information. 

Note that by adding together Eqs. (2.21) and (2.22) all of the odd-spin terms cancel 
and we recover the crossing symmetry constraint for the operators appearing in the Re[0] x 
Re[0] OPE. Alternatively, we could subtract Eq. (2.22) from Eq. (2.21) to obtain a crossing 
constraint which relates just the odd-spin operators appearing in x 0* to the operators 
appearing in x 0. Also notice that the (— I) 1 factor in Eq. (2.22) cancels against the (— 1)' 
factor that occurs in the definition of the conformal blocks, so that in this equation the odd- 
spin terms are qualitatively similar to the even-spin terms. This may be contrasted with 
Eq. (2.21), where odd-spin terms have the opposite sign relative to the even-spin terms. For 
this reason, we have found that it is much easier to obtain a well-behaved linear program 
using the constraints of Eq. (2.22) as compared to the constraints of Eq. (2.21). Thus, in 
the present work we will mainly focus on the bounds that can be obtained using Eq. (2.22), 
though in future studies it may be useful to incorporate the full set of constraints. 

Let us also briefly mention another way to generalize the procedure outlined in the 
previous section. Thus far, we have only used the knowledge that the unit operator appears 
in the crossing relation, but in many situations one might have additional information. For 
example, if it is known that an operator O of dimension A and spin / appears in the 0x0 
OPE in addition to the unit operator, and we also know its three-point function coefficient 
Xq, then one can simply make the replacement 1 — > 1 — X 2 ~F^j in the objective function 
of the linear program. This modification can then lead to more stringent bounds. It is 
particularly straightforward to implement in the case of the stress tensor T ah or a conserved 
global symmetry current J a , since in these cases the dimensions are known and the A's are 
fixed by Ward identities, as we will review in Section 4. 
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3 Superconformal Blocks 



At this stage we could proceed to derive bounds on 3-point function coefficients in conformal 
field theories with global U(l) symmetries. However, because we would also like to derive 
similar bounds in M = 1 superconformal theories, we will first consider more carefully the 
additional constraints imposed by supersymmetry. In particular, three-point functions of 
primary operators in the same supersymmetry multiplet are related to each other by the 
superconformal algebra, and one can construct "superconformal blocks" which sum up the 
contributions of all operators in a given superconformal multiplet. 

We will focus on four-point functions involving a complex scalar that is the lowest 
component of a chiral superfield $ of dimension d = §-R<j>. In terms of the operators appearing 
in the x 0* OPE, the superconformal block decomposition looks like 

(0(Z1)0*(Z 2 )0(Z3)0*(Z4)> = l A o| 2 (-l)W^)- (3-1) 

12 34 0g $ x $t 

Here, we have adopted the notation G $ x $t to indicate that the sum is over super- 
conformal primaries O appearing in x 0*, and not simply primaries under the conformal 
subgroup. By definition, superconformal primary operators O are annihilated by the S 
and S generators in the superconformal algebra, from which it follows that they are also 
annihilated by the K generator. However, a finite number of superconformal descendants 
of O are also killed by K, so one may decompose Qa,i(u,v) into a finite sum of conformal 
blocks gA,i(u, v). 

Just as the explicit expression (2.8) for conformal blocks was crucial for the analysis of 
[33-35], an explicit expression for superconformal blocks will be crucial for us. We find that 
J\f = 1 superconformal blocks in the 0x0* channel are given by 

_ (A+j) (A - / - 2) 

Ga,i - 9a,i ~ 2{A + l + 1} SA + i,m - g(A _ l _ 1)9 a + i,i-i 

(A + /)(A-/-2) 
+ 16(A + / + l)(A-/-l) fi ' A+2 ' z - (3 ' 2) 

To our knowledge, this expression has not yet appeared in the literature, though analogous 
results for M = 2 and M = 4 theories are known [37]. Eq. (3.2) is the key ingredient we 
need to apply the technology of Section 2.2 to superconformal theories. In the following 
subsections, we will give two derivations — one involving explicit analysis of superconformal 
two- and three-point functions, and another quicker but less illuminating argument leveraging 
known expressions from M = 2 theories [37]. The discussion is somewhat technical, and 
readers interested solely in bounds on dimensions and OPE coefficients should feel free to 
skip to Section 4. 

Our first derivation of Eq. (3.2) proceeds as follows. We start by understanding which 
superconformal primary operators O ai '" ai can appear in the OPE 0x0*. We then determine 
which superconformal descendants of O ai " ai are conformal primaries, and further calculate 
the relationships between two- and three- point functions of these conformal primaries. Since 
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each conformal primary contributes a block gA',i> to (00*00*), we can piece together Q^,i 
from these contributions. For completeness we also include a brief discussion of the 0x0 
channel. However, in this case only a single operator in each supersymmetry multiplet 
may contribute, so the superconformal blocks turn out to be the same as the conformal 
blocks Eq. (2.8). Our conventions for the superconformal algebra and spinor notation are 
summarized in Appendix A. 

3.1 Superconformal Three-Point Functions 
3.1.1 0x0 OPE 

Let us start by examining the 0x0 OPE, since the constraints from superconformal symmetry 
are particularly transparent in this case. This analysis is not needed later, but we include it 
for completeness and to establish some notation. For some previous discussions of this OPE, 
see [46, 47]. In this subsection we will follow the notation and conventions of [45], where 
a superconformal primary O 1 (I denotes Lorentz indices) is specified by spins and 
conformal weights (<?e>, <?£>), w hich are related to the dimension and i?-charge via qo + ~q~o = 
Aq and \{qo ~ Qo) = Ro- The unitarity bound for non-chiral superconformal primary 
operators then requires [42-44] 

A / > I^V-J+Jl+J+J + 2. (3.3) 

To begin, note that since Q<f){x) = 0, only operators that are annihilated by Q may 
appear in x 0. A priori, there are four possibilities: 

1. Chiral primaries. Since these transform in (j,0) representations of the Lorentz group 
SU(2) x SU(2), they can appear only if j = 0. We will denote the linear combination 
of chiral primaries appearing in x by 2 . 

2. Descendants of the form Q ai Q°®~ a i) a i- a i ^ where / is even and O 1 satisfies the short- 
ening condition Q^Q^-m^-^i = q. (Note that this implies (fo 1 = 0, so that QO 1 
is indeed killed by Q.) The superconformal algebra implies [45] that such operators 
satisfy q — (I + l)/2. Then using Rq i = 2i?$ we find A i — 2d + I — |, so that the 
dimensions of these operators are determined by their spins. We will denote the linear 
combination of these descendants with spin / as QO\ . 6 

3. Descendants of the form Q & Q aa i- a i a i- a i^ where O satisfies the shortening condition 

Qa 2 ...ai +2 )a 1 ...ai _ g_ g^jj multiplets must satisfy q = —(I + l)/2, which implies 
upon matching i?-charges that = 2d — I — 5/2. However, this violates the unitarity 
bound Eq. (3.3), so such operators actually cannot appear. 

4. Descendants of the form Q Q 2 ^ n O I , with n = 0, 1, 2. 

6 We are grateful to Alessandro Vichi for pointing out the possibility of these operators in the 0x0 OPE. 
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Thus, we expect the OPE to take the form 

mm = C(x,P)0 2 (O)+ J2 CifaPffiOlW + ^tfCfaP^O'tp), (3.4) 

1=2,4,... O 1 

where the latter sum runs over superconformal primaries with Rqi = 2/2$ — n, and a priori 
n — 0, 1, 2 depending on how many powers of Q appear. 

We can obtain additional constraints on the operators O 1 by acting on both sides of 
Eq. (3.4) with an S generator. Note that S kills the left-hand side because [S, P] ~ Q and 

is chiral and primary. On the right-hand side, we can commute S through all powers of 
Q and P, since {S, Q} = and Q [S, P] ~ Q = 0. However, if powers of Q were present, 
there would be terms involving {S, Q} which would not vanish when acting on O 1 . Thus, 
we conclude that Ci(x, P,Q) = Cj(x,P) and therefore Rqi = 2i?$ — 2. In this case the 

1 indices must correspond to even-spin operators due to the symmetry under exchanging 
x -H- —x. Finally, the unitarity bound Eq. (3.3) implies ^q 2 i > |3i?$ — 3| + 1 + 3. Note also 

that Q O 1 is primary under the conformal sub-algebra. 

Instead of playing directly with the superconformal generators, an alternative approach 
that will prove useful later is to consider the general form of superconformal-covariant three- 
point functions. Let us take a moment to recover the above results using this language. 

The $ x $ OPE contains a superconformal multiplet O 1 if and only if the three-point 
function {&(zi+)§(z2+)0 1 ^ (z^)) is non-vanishing, where the z's are superspace coordinates 
(x, 9, 9), and z + indicates dependence only on the chiral subspace (x + i9a9, 9). The general 
form of such a three-point function consistent with superconformal symmetry is 

(4>( 2l+ )4(, 2+ )0'» W > = t,{X l^ Z \ (3.5) 

X 31 32 

where Xjj = £j_ + 2i9j<j9i — Xj + denotes the supertranslation-invariant interval built out of 
anti-/chiral coordinates Xi± — Xi ± i0i<T9i, X 3 and ©3 are given by 

"I oc^ oc^ 

*s = -- % T2 2 2 H r(a a a b a c a d ), (3.6) 

1 X \3 % 32 

©3 = i^f(r a 9 31 -i^f-a a 9 32 , © 3 = 6j, (3.7) 

X 13 X 23 

and t 1 has the homogeneity properties 

^(AAX 3 ,A©3,A© 3 ) = A 2a A 2 V(X 3 ,03,©3) (3.8) 
with a = \{2qo + q — Ad) and a = \{qo + 2g G — 2d). 

Since the covariant derivative vanishes when acting on the left hand side of Eq. (3.5), 
we obtain an additional constraint (using Eqs. (6.1) and (6.2) in [45]) 

= ^ / (X 3 ,0 3 ,©3) 
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which implies that t 7 (X 3 , 63, 63) = t 1 (X 3 , 63), where X 3 = X 3 + 2i0 3 a0 3 . Finally, under 
Z\ -H- z 2 we have X 3 <H> — X3 and 83 -h- —63. There are three possible solutions to these 
constraints, 

t I (X 3 ,e 3 ) = const., (3.10) 
corresponding to O 1 being a chiral "$ 2 " operator with Rq = 2/?$, 

^(x 3 ,e 3 ) cx ef 1 xri 2 ...xri = ®f x x 3 %...x 3 % (3.11) 

corresponding to the short operators Of, and 

t T (x 3 ,e 3 ) oc egXf " 2 ^" 1 ^ 1 ...^ = Qlx£ o - 2d - l - 1 x%\..xi i , (3.12) 

corresponding to O 1 being a non-chiral operator with Rq = 2/2$ — 2. Since the only 
irreducible Lorentz representations that can be built out of a single vector X 3 (or X$) are 
traceless symmetric tensors, O 1 = O ai '" ai must have definite integer spin I = 2j = 2j, and 
invariance under z\ -H- zi further tells us that I must be even. The descendant operator Q O 1 
then has the correct quantum numbers to appear in the 0x0 OPE, in precise agreement 
with the preceding argument. 

Here we see that for each supermultiplet appearing in $ x $, there is exactly one 
conformal primary appearing in x 0. This is essentially because 2 , QOf , and Q O 1 are 
the only conformal primaries in their respective supermultiplets with the correct i?-charge. 
Consequently, the superconformal blocks for decomposing (00*00*) in the 0x0 channel are 
the same as the conformal blocks. Next we will turn to considering the 0x0* channel, where 
this will no longer be the case. 



3.1.2 0x0* OPE 

We determine which operators can appear in the 0x0* OPE by examining three-point 
functions (^^O 1 ^). Once again, let O 1 be a superconformal primary with conformal weights 
(qo,Qo) an d spins Following [45], we must have 

($(z 1+ )$t(^)o 7 t(^)) oc _l_^(X3,e 3 ,e 3 ), (3.i3) 

2; 31 X 23 

where t 1 satisfies Eq. (3.8) with a = |(2go + q ) — d and a = \{2q + q ) — d. 

Demanding the appropriate chirality properties imposes further constraints. Just as in 
the 0x0 case, requiring D^t 1 = means t 1 must be a function of X 3 and O3. We must 
additionally require 

x 61 - 8 

= Dtf{X 3 ,Q 3 ) = i^—t\X 3 ,® 3 ), (3.14) 

x 23 oe 3 

so that t 1 is actually a function of X 3 alone. Note that since the i?-charge of X 3 vanishes, 
the i?-charge of the correlator ($$^(9 7 ^) must vanish as well, which means O 1 = O i>[ should 
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be a real operator with qo =~q~o- Since we again can only build Lorentz representations out 
of a single vector X 3 , the only possibilities are traceless symmetric tensors, so O 1 = O ai '" ai 
must have definite integer spin I — 2j — 2j. 

In summary, we have found that the only superconformal primaries appearing in the 
$x $t OPE are traceless symmetric tensors O ai '" ai with vanishing i?-charge. Superconformal 
symmetry determines the 3-point function to be 

mz 1+ )&(z 2 „)(D"i(z 3 )) oc -^X^'^Xl 1 ...^ -traces. (3.15) 

In this case, the unitarity bound Eq. (3.3) requires Aq > / + 2. 7 The operators which can 
enter the OPE of the lowest components <px <fr* are then i?-charge zero descendants of a real 
superconformal primary, P n (QQ) m O ai '" ai . To understand how these operators contribute 
to the four-point function (</></>* 4><f>*), we must now organize them into representations of the 
conformal sub-algebra. 



3.2 Decomposition of Superconformal Multiplets into Conformal 
Multiplets 

In this section, we will examine the structure of a multiplet built from a real superconformal 
primary O ai " a ' of dimension A. The full superconformal multiplet can be decomposed into a 
direct sum of conformal multiplets, connected together by supersymmetry transformations. 
Here we will show explicitly how this decomposition works for operators that appear in the 
(f) x cf)* OPE — namely operators of vanishing i?-charge and definite spin. As a result, we 
will see how superconformal symmetry relates the OPE coefficients of different conformal 
primaries, and consequently how Qa,i decomposes into a sum of <7a/s. 

Note that O ai '" ai is symmetric and traceless in its indices. Throughout this subsection, 
we will adopt the convention of implicitly symmetrizing and subtracting traces in a,i for 
i = 1, . . . , I. This has the virtue of greatly simplifying notation, though one must be careful 
when manipulating expressions. 

A convenient way to describe the descendants of a superconformal primary operator 
O ai '" ai (0) is through superspace. For example, defining the superfield O ai '" a '(x, 8, 8) = 
e xP+BQ+eQQa 1 ...a l ^) ) ^ we j^yg ^g component expansion 

O ai - ai {x,8,8) = A ai - ai {x) + ( a B aai - a '{x) + ( 2 D ai - at {x) + ... (3.16) 

where ( a = 8a a 8, and ". . . " represents fields with non-zero i?-charges. The component fields 
■Qaa x ...ai an( j ^ai...a, are ^j ien re i a t e d to A ai '" a; through the action of Q and Q as 

j^aa\...ai _ _^^a j^.\...ai yj^ 

D" 1 = -—E a B aai - ai - —d 2 A ai ~ ai , (3.18) 
64 16 V ' 

7 With an exception, of course, for the unit operator which has A = I = 0. 
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where we have defined S c 



[Qai Qc 



Both A ai '" ai and D ai ' ai are in the spin-/ representation of the Lorentz group, but B aai,,,ai 
can be further decomposed into irreducible representations. Recall that under SO (4) = 
SU(2) x SU(2), the spin-/ representation of SO(4) transforms as (j, j) with j = 1/2. Since 
B aCLl --- a ' has an additional vector index, it transforms as 



;i/2,i/2)®(j,j) 



(j + 1/2,j + 1/2)©(j-1/2,j-1/2) 
© (j + 1/2, j - 1/2) © (j - 1/2, j + 1/2) . 



(3.19) 



The first two components on the right-hand side are a spin-(/ + 1) representation J aa ^--- a i = 
^(aai...a ; ) _ t raceS; anc [ a spin-(/ — 1) representation jV a2 '" ai = B ' b bai . The remaining two 
components comprise an operator £ a °i-°a which is traceless and has vanishing total sym- 
metrization. L can be further decomposed into irreducibles by projecting onto its "anti/self- 
dual" parts, satisfying £ J '^ 1 - ai — +ij^e aai bc L b ± a2 ''' ai (although this will not be important in 
our discussion). Notice that since L is not in a traceless symmetric representation, a primary 
operator built from it cannot appear in the OPE of <fi with <ft*. Nonetheless, it will play a 
role in the identification of conformal primaries below. Altogether, we may write 



B 



J 



aa\ ...cli 



+ 



l 2 



(I 



(12... a; 



+ L 



aa± ...0,1 



(3.20) 



where as usual we are implicitly symmetrizing and subtracting traces in the a^. The coeffi- 
cient of N is such that the projection N a2 " Ml = B b bai " ,ai works correctly. 

Now let us consider the action of a special conformal generator K a on the components 
of O. We will be interested in determining which linear combinations of superconformal 
descendants are annihilated by K a . After some algebra, one can determine the action 



K„ 



j^ba!...ai 
^bai ^pcj^da2...a { 
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2(A-1) 



(e ba \ d A da *- a <) 



(3.21) 



as well as 



Kn 



( 



Jjai...a t 
p2j^a\...ai 
p pat j^ba 2 ...a l 



\ e aib cd P b B cda i- a ' J 



4(A 

C 

V 2(J- 




-4/ 
2(A-/-2) 
-2(1 - 1) 





4/ 

2(A + l) 
-2(1 + 1) 



2 





2(A-1) J 



( P a A a i- a ' \ 
P ai A a a2 "' ai 
5Z 1 P b A ba *- ai 
K e a \ cd B cda *- a > J 

(3.22) 



from which we find that the linear combinations 

/ 



-naa 1 ...ai 
prim 



T-jOi...a; 

prim 



paa\...ai 

p,ai...ai 
I 



A - 1 
1(1 + 1) 



^aaj ^pc j^da 2 . . .a t 



'A(A-l) 



Aai—ai 

8(A- l) 2 

e aib cd P b B cda2 - ai 



4(a-i; 



_p pai j^ba2-..ai 



(3.23) 



(3.24) 
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are primary operators under the conformal subgroup. Note that only the L component of B 
is shifted in the above expression for -B prim , so that J and N are already primary. 

An important fact is that when the unitarity bound A > / + 2 is saturated, our super- 
conformal multiplet is "shortened," and the descendants N, L prim , and -D pr i m actually vanish. 
For example, the supercurrent J a (z) with A = 3 and I = 1 contains only the i?-symmetry 
current J%{x) and stress tensor T ab (x) as conformal primary components with vanishing 
i?-charge. This will be reflected in explicit calculations below. 



3.3 Conformal Primary Three-Point Functions 

Next we would like to see how the three point functions (<jxf)*J), (<fxj)*N), and D pT - im ) 
are related to ((jxj)*A). We will also verify that (00*L prim ) = 0, as expected because L prim is 
not in an integer-spin (traceless symmetric) representation of the Lorentz group. 

Let us set 9\ = 8 2 = Q\ = 8 2 = 0, and 8 3 = 8,8 3 = 8 in the correlator Eq. (3.15) to get 
the 3-point function (<fr(xi)<j)*(x 2 )O ai "' a ' (x 3 , 9, 9)). Next, expanding in 9,9 and comparing 
with our component expansion Eq. (3.16), we find 

*A ai - ai ) = A-z A-i Zai --- za ' ( 3 - 25 ) 

x 13 x 23 

A-2q-l 

*ja ai ...a^ = j(A + £) 12 1 A _ t Z a Z ai ...Z ai (3.26) 

x 13 X 

A-2q-l 

' L aai...at^ = 2/ ^ A _, Y aai Z a2 . . . Z a > (3.28) 

x 13 X 23 

A-2q-l /, ai 1(1-1) 

*D ai - ai ) = A ^ A _, -Z 2 ^-Z a2 ...Z ai + ^ J - Vbc Y aib Y a2C Z a3 . . . Z ai (3.29) 

X 13 X 23 V 2 x 21 ^ 



x 13 X 23 

(^*N a *- a ') = / A ~ Z ~ 2) ^ + 1) ^ A _ l Z 2 Z a * . . . Z a > (3.27) 

M x 13 x 23 



'13 23 

'(A + 0*32-3:13 , (A + /)(A-/-2) 



x 13 x 32 § 



where 



31 Xr, 



ya _ ^31 _ A 32 v ab — r ab ,r c r d (1 

Zj — „ „, I — 22 t cdX 3 \X 32 , yo.ov) 

,-y» ij ry^ ry^ ry^ 

■ L 31 x 32 J '32 J '31 

and we are implicitly projecting the right-hand side of each expression onto the appropriate 
Lorentz representation (symmetrizing and subtracting traces as necessary). Using Z 2 = 
^12/(^31*32)) we see that the correlators (4><f)* A), (4>(fi*J) and (00* iV) take the expected form 
for a 3-point function of conformal primary operators. Further, taking the appropriate 
derivatives of the above expressions and constructing the linear combinations corresponding 
to L prim and -D prim , we obtain 

L^ a ') = (3.31) 
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as expected, and 

Urh* D ai - ai \ A(A + /)(A - / - 2) X 12 q 7 2yax ?ai fo oo\ 

- x ^ x ^ zz ■ Z ■ (3 - 32) 

Notice that three-point functions involving N and -D prim vanish when A = 1 + 2, which is 
precisely what we expect for short multiplets that saturate the unitarity bound. 



3.4 Conformal Primary Norms 

Finally we must determine the normalization of the two-point functions (J J), (NN), and 
(-Dprim-Dprim) • One could do this either by expanding out the superconformally covariant 
expression for the two-point function of O derived in [45] into its various components, or by 
using the explicit expressions for J, N, and -D pr i m in terms of Q, Q, and P acting on A, and 
using the superconformal algebra to compute their norms in radial quantization. We here 
adopt the latter approach. We refer the reader to [44] for many examples of this type of 
computation. 

To begin, we assume that the superconformal primary operator A is canonically normal- 
ized 

(A bl - bl \A ai - ai ) = symmetrize^ 1 61 . . . rj a ' b ') - traces 
= - V aiK(1) ■ ■ ■ V aiK(l) ~ traces 

_ 2 * 1 .„a, i 6i...6, j (333) 

where we've defined X^ 1 '" ah 1 "' 1 for future convenience, and \A ai - ai ) = A ai "' a '(0)|0) is the 
state created by the operator A ai " Ml (0) in radial quantization. 

Next we would like to determine the normalization of B^^' ai . Starting from Eqs. (3.17) 
and (3.23) and working through the algebra, we find that 



(B^ bl \B^) = 2((A(A+l)-/ 2 -^^)r / X 1 +^(2A + 2/ + l + ^ r )^- 

-/ \2A-2l + l- ^— yj n h cl V aai j (A bl - bl \A Cia2 - a >), (3. 
from which we can extract the component normalizations 

(jbh-b^jaaL.-atj = 2(A + /)(A + / + l)2^""" a!;66l '" 6i , (3.35) 



as well as 
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where we have used the relation r j ah x^ a2 --- ai ' bb2 --- bl — (i+l)_x^' 1 ' ai ' b2 '' Ml . Although we will not 
need it, for completeness we also have 

I T bb 1 ...b l \ T aa 1 ...a l \ _ 8/ 2 A(A + /)(A — I — 2) a 6 T ai...q,;6i...6; /o q 7 \ 

Wprim l-^prim / — (/ + 1) 2 (A — 1) ' \O.Ol) 

where we are implicitly subtracting traces and the full symmetrization (in either the a, a; or 
b, hi indices) from the right hand side — that is, projecting onto the Lorentz representation 
corresponding to L. 

Finally we must determine the normalization of -DpJ im az • In order to simplify the calcula- 
tion, it will be helpful to write everything in terms of primary fields, 

n a i--- a i _J_^ r>aai-(H f(j + jj + ~ jj + jj p2 401. ..a; , ^ pai p 4602-0; 

prim 64 ^P"m 16(A-1) 2 + 8(A-1) 2 6 

q; 

- 16(A _ 1) ^ l6 c^S- ai (3-38) 

so that 

64 2 (^-''|B°; lm '") = (BX^m^B^-*) - (^^{B^'KSO'P-'P^ -">) 

+4 i(i + 1) ^_l )(A + 1) ( ^ |(5t)tp2| ^ a , ) 

12/ 66 6 d 

+ (A - l) 6 1 cd ^P r ' m 'K" 6 ) -^el-^prim /) (3.39) 

where we have used that all terms of the form ((. . . )K\D pr - im ) vanish. Evaluating each of 
these terms using the superconformal algebra and putting everything together, we obtain 
the final result 



mb ,.. blln a,.. ai) A 2 (A - / - 2)(A - / - 1)(A + Q(A + I + l)_ ai ... aiM ... bl 

X-^prim l-^prim / — _ I ' 



3.5 J\f = 1 Superconformal Blocks 

To summarize the results in the previous subsections, we have found the three-point function 
coefficients 

= 1 

= i(A + 

,(A-/-2)(/ + l) 





<t>*A 


\ 


></>* J 






A<A, 





2/ 

A(A + 0(A-Z-2) 



8(A-1) (3 ' 41 > 
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and the norms 

(A\A) ~ 1 

(J|J) ~ 2(A + /)(A + / + 1) 

(iV|iV> ~ 2(Z + im-Z-2)(A-f-l) 

- ^-'-^-_^ + 0(A + / + l) f (342) 

where "~" means multiplied by the appropriate canonically normalized tensor. Combining 
these results, we find the dimension A, spin I superconformal block given in Eq. (3.2), which 
we reproduce here for the reader's convenience, 

_ (A + Q (A-/-2) 

(A + Q(A-/-2) 
+ 16(A + / + 1)(A-/-1)^ A+2 ' Z - (3 ' 43) 

A few comments are in order. First, I = is special, since in this case the N component 
does not exist. However, one can consistently take <7a,-i = 0, and then the above equation 
correctly accounts for this situation. Second, in the case of superconformal primary operators 
that saturate the unitarity bound, A = 1 + 2, the third and fourth terms vanish, which is 
precisely what we expect due to the fact that the N and D prim components are not present 
in short multiplets. Finally, in the case of the unit operator, with A = I = 0, the second and 
fourth terms vanish due to the coefficient going to zero, and the third term vanishes because 
the conformal block goes to zero. Thus, we simply obtain that £/ ,o = 9o,o — 1- 

Let us also note that Eq. (3.43) determines the superconformal blocks for four-point 
functions of all component fields in not just the lowest component <f)(x). The reason 

is that there are unique superconformally-invariant extensions of the conformally-invariant 
cross-ratios u, v with the correct chirality properties to appear in a four-point function 

($(^i + )$ t (22-)$(2 : 3+)$ t (-24-))- They are given by [45] 

« = :Jh?> ^ = V41 1 )' ( 3 - 44 ) 

X 23 X 41 Z 



where the x's in the trace should be thought of as bispinors, [x) aa = x a a c ^ x and {x 
—XaO^jx 1 . Since u and v become u and v when we set Q{ = 6i = 0, we must have 



act 



($(^ 1+ )$t(^)$(^ + )<f t (^4-)) = X lG±Mv), (3-45) 

12 34 c» g $ x $t 

where ^a,z is given by Eq. (3.43) above. One can now perform 9, 9 expansions on both sides 
to derive the superconformal blocks for specific component fields. 

Finally, let us mention that it may be possible to derive the superconformal blocks by 
mimicking the derivation of g^i in [41]. One would start with the expansion Eq. (3.45) and 
apply the quadratic casimir of the superconformal group acting on $(zi + ) and $(22-) to 
obtain a differential equation for Qa,i, which could then be solved. 
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3.6 Deriving M = 1 Blocks From M = 2 Blocks 



In [37], Dolan and Osborn computed superconformal blocks for four-point functions of a 
particular kind of BPS operator in M = 2 theories, using Ward identities special to higher 
supersymmetry. At the very least, we should be able to decompose their expression into 
M = 1 superconformal blocks Qa,i- However, requiring that this is possible gives a strong 
consistency condition on Q A \ — so strong in fact that it determines Q A) i completely! In this 
subsection, we will use this fact to give an alternate derivation of Eq. (3.43) that requires far 
less computation than in Sections 3.1-3.5, though it leverages important results from [37]. 

The operator ip 1 ^ considered in [37] is a triplet under SU(2)#, neutral under U(1)r, and 
has scaling dimension 2 (here i,j = 1,2 are SU(2)# indices). It satisfies the BPS conditions 
Q£tpi k ) = e^Q^if^^ = 0, which imply that under the H — 1 sub-algebra generated by 
and Qa\, the operators ip u ,ip 21 , and p 22 are anti-chiral, linear, and chiral respectively. The 
important fact for us is that ip 22 = <fi is chiral, so ((fxp* <$<$*) can be decomposed into a sum 
of Ga/s- Note that the form of Q A) \ is independent of the dimension of <fi. In particular, it 
is irrelevant for our purposes that <fi is restricted to have dimension 2. 

Any M = 2 multiplet that can appear in the OPE p 1 ^ x p kl must be built from a primary 
of dimension A and definite integer spin /. We will denote such a multiplet by (A)^ =2 . 
The "extra" supersymmetry generators Q 2 , Q 2 connect different M = 1 multiplets within 
(A)^ =2 exactly analogously to the way Q and Q connect different conformal multiplets 
within (A)^ =1 , as discussed in Section 3.2. Thus, we have the decompositions 

(A)^= 2 = (A)f =1 © (A + l)£f © (A + 2)f =1 (3.46) 
(A)^ =1 = (A)f= ©(A + l)^T°©(A + 2)f= , (3.47) 

where we have ignored multiplets which cannot appear in the OPE of two scalars. We can 
then write the ansatze 

<F 2 = gA,i + N{A,l)g A+1>l ^ + J{A,l)g A+1>l+1 + 0(^^+2,1 (3.48) 
Ga,i = 9a,i + n(A, Z)g-A+i,«-i + j(A,Z)#A+i,m + d(A,l)g A+2 ,i, (3.49) 

where N, J,D,n,j,d are functions we would like to determine. Note that j, n, and d must 
be rational functions of A and I. This is clear without any computation, simply from the 
viability of our first method for determining Q A \ (Sections 3.1-3.5). 

Using formulae from [37], we find that the M = 2 superconformal block contributing to 
00*) is given in terms of conformal blocks by 

VA,l = 9A,l - gA+1,1+1 - ~9A+l,l-l + jQA+2,1 

(A + / + 2) 2 (A + / + 2) 2 

4(A + / + 1)(A + I + s)^ 2 ^ 2 16(A + / + 1)(A + I + 3)^ A+3 < Z+1 

(A-/) 2 (A - If 

+ 64(A - I - 1)(A - I + i) 9A+2 > l ~ 2 64(A - I - 1)(A - I + i)^ 3 -'" 1 

, (A + / + 2) 2 (A-/) 2 

256(A + I + 1)(A + I + 3) (A - I - 1)(A - I + i) 9A+4 ' 1 ' {6 - 0U) 
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Upon comparison with Eqs. (3.48) and (3.49), each coefficient in the above expression 
implies an equation relating N, J, D,n,j, and d. We will solve these equations by first 
determining j and n, and finally computing d in terms of them. To begin, the <7a+i,j+i and 
9a+2,i+2 terms in Eq. (3.50) imply 

-1 = J(A,/) +J (A,0, and ^^rj^f^ = J(A, Z)j(A + 1, Z + 1). (3.51) 
With some foresight, but without loss of generality, let us make the substitution 

J'(A,0 = ~ 2 (A + I + 1) (1 + " (A + I, A - Q), (3.52) 
where a(x,x) is a rational function we must determine. Then Eqs. (3.51) imply the equation 

/ — \ x + 2 a{x + 2,x) 

a(x,x) = ; r, (3.53) 

V ' ; x l + a{x + 2,x) K J 

and it's not difficult to show that any rational solution a(x,x) must vanish identically. 
Consequently, we obtain 

A similar analysis using the gA+i,j-i and gA+2,1-2 terms in Eq. (3.50) gives 

n (A, Z) = - (A ~ / ~ 2) N , N(A, I) = - ,[ A ~ l) - . (3.55) 
v ' ; 8(A-/-1)' v ' 1 8(A-Z-1) V ; 

Finally, let us solve for d(A, I). The gA+4,1 term in Eq. (3.50) determines D(A, I) in terms 
of d(A + 2,/). Plugging this in, along with our solutions for N, J,n, and j, the remaining 
terms in Eq. (3.50) imply equations with the following structure 



9A+2,l 
<?A+3,i+l 



d{A,l) ~ d(A + 2,Z) (3.56) 
d(A + 1,/ + 1) ~ rf(A + 2,Z) (3.57) 
d(A+l,Z-l)~d(A + 2,Z), (3.58) 



where "~" means "is algebraically related to." Making the substitutions A — > A — 1 and 
Z — > I — 1 in Eq. (3.57), we are left with three algebraic equations relating three "variables" 
cZ(A, Z), d(A + 2,1), and cZ(A + 1, Z - 1). Solving them gives 

(A + Z)(A-Z-2) (A + Z + 2)(A-Z) 

rf(A ' /)= 16(A + Z + l)(A-Z-ir D(A ' /) = 16(A + Z + 1)(A-Z-1) - (3 ' 59) 



To summarize, we have re-derived Eq. (3. 43), 8 and also obtained the decomposition of 

8 It's possible that similar arguments suffice to determine M = 2 conformal blocks from N = 4 conformal 
blocks. If this is the case, it's fascinating that a maximally supcrsymmetric result, which can be derived 
using special properties of M = 4 BPS multiplets, completely determines the corresponding results for lower 
supersymmetry. 
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M = 2 conformal blocks into M = 1 conformal blocks 

M=2 _ (A + l + 2) (A-Q 

~ ^ " 2( A + / + 1) " 8( A - / ^ ^ A+1 

(A + / + 2)(A-/) 



i6(A + / + i)(A-/- r 



^A + 2,i- (3.60) 



4 Bounds 



Now we finally turn to using the results obtained in Sections 2 and 3 to obtain bounds on 
CFT and SCFT data. We will start by considering bounds on the OPE coefficient of the 
lowest-dimension scalar appearing in the $ x $t OPE (which we call "$^$"), where $ is a 
chiral multiplet in an M = 1 superconformal theory. When the dimension of <3> is somewhat 
close to 1, we find that these OPE coefficient bounds are sufficiently strong to yield an upper 
bound on the dimension of This is a completely general result about the dimensions of 
non-chiral operators in strongly-coupled M = 1 superconformal field theories. We will also 
present a bound on the OPE coefficient of an arbitrary scalar operator that can appear in 
this OPE, independent of any assumptions about the spectrum. 

Then we turn to bounding the OPE coefficients of flavor currents. In general CFTs these 
are spin-1 operators J aI of dimension 3, and in M = 1 theories the J aI are embedded into 
real scalar operators J 1 of dimension 2. We will review how Ward identities fix these OPE 
coefficients in terms of the coefficient of (J 1 J J ) oc r IJ and the charges of 0, allowing us to 
bound the quantity tijT i T j , where r/j = (t /j ) _1 and T 1 are the generators of the flavor 
symmetry in the representation. Roughly speaking, t ij measures the number of degrees 
of freedom charged under the global symmetries, 9 and our bound says that the effective 
number of degrees of freedom that are charged cannot be much smaller than 1. We present 
this bound in both non-supersymmetric and supersymmetric CFTs. 

Finally we consider the OPE coefficient of the stress tensor, which is similarly fixed by 
Ward identities in terms of the dimension d of <fi and the central charge c. This will allow 
us to derive a lower bound on the value of the central charge in both non-supersymmetric 
and supersymmetric CFTs. In the former case, the stress tensor is a spin-2 operator of 
dimension 4, and the bound will assume only the existence of a real scalar primary operator 
of dimension d. In the latter case, the stress-tensor is the 6a a 6 component of a spin-1 U(1)r 
current multiplet of dimension 3, and the bound will assume only the existence of a chiral 
primary scalar of dimension d. 



4.1 Dimension of and Scalar OPE Coefficients 

We start from the crossing relation Eq. (2.22), which involves only the <fi x <ft* channel of 
the four-point function (0^*00*). Superconformal symmetry additionally allows us to group 

9 For example, if the flavor symmetry is weakly gauged with coupling g, then t ij is proportional to the 
contribution of the CFT to the beta function of g. 
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terms into superconformal blocks, so that we may write 

£ |A o | 2 (-l)Wu,t0 = E l^| 2 (-l)^A,/(t',n), (4.1) 

Oe*x$t Oe3>x$t 

where as before we have written $ x $^ instead of x 0* to indicate that the sum is over 
superconformal primaries in the OPE of <fi with </>*, as opposed to simply primaries under 
the conformal subgroup. 10 

From here, the procedure is exactly as described in Section 2.2. Suppose the operator 
Oq with dimension Ao is the lowest-dimension scalar appearing in $ x Isolating the 
contributions of Oq and the unit operator, we have 

|Ao | 2 .Fao,Jo = 1 - /2 \^o\ 2 Fa,u (4-2) 

where J-a,z is given by Eq. (2.11) with — > ( — Finally, to obtain the best possible 

bound | Ac? 1 2 ^ a (l); we must minimize a(l) over all a C V* satisfying the constraints 

• o;(Ja,o) > for all A > A , 

• a (-Fa,/) > for all A > / + 2 and I > 1 (not necessarily even), 

• a(JAo,o) = 1- 

Then if the resulting bound tells us that |Ae> | 2 < a(l) < 0, there is a contradiction with 
unitarity, and we learn that Oq cannot have dimension A . 

Let us highlight some assumptions implicit in this procedure. Firstly, we are assuming 
that $ is uncharged under any global flavor symmetries (that is, non-i? symmetries), since 
otherwise there would be a symmetry current J of dimension 2 in the OPE $ x $t, which 
would necessarily be the lowest-dimension scalar by the SUSY unitarity bound. Alterna- 
tively, if $ has flavor charges, and we wish to bound the lowest-dimension scalar in $ x $^ 
that is not a flavor current, then we must incorporate flavor current blocks J-2,0 m to the 
objective function of our linear program 1 — > 1 — lAjpJ^.cb as discussed in Section 2.3. 

Secondly, note that we are only using part of the full crossing relation Eq. (2.20), and it 
is possible that one could obtain stronger bounds by incorporating the additional relations 
Eq. (2.21) (whose terms also can be grouped into superconformal blocks if desired). So far, we 
have not had success incorporating these extra constraints into a well-behaved linear program 
- namely one where our choices of finite-dimensional subspaces V\4 and discretizations 
D = {(Aj, li)} lead to answers that don't violate other constraints Fa>.i> > for (A', /') ^ D. 
It is certainly possible that these difficulties can be circumvented. However, in this paper, 
we choose to focus on the information that can be learned from Eq. (4.1). 

10 Note that the methods in this section apply equally well to a general CFT with a global U(l) symmetry 
in which case the bounds are strictly weaker. 
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Figure 4 shows the resulting bound on the dimension of <3>^<3> as a function of the dimension 
d of <3?. Here we have taken k — 6, and then for each value of d we scan over values of A 
(with a spacing of 0.01) until we find the smallest dimension such that there is a contradiction 
with unitarity. As d — > 1, we see that the bound approaches 2 from above, consistent with 
the existence of the free theory. (Bounds very near d — 1 are computationally intensive 
to obtain, so we defer very close exploration of this region to future work.) On the other 
hand, we see that the bound becomes very weak and shoots off to infinity around ~ 1.16. 
For dimensions larger than this value, the resulting bounds on |Ae> | 2 become stronger and 
stronger as A becomes large, but never lead to a violation of unitarity. We also note that 
at k < 6 we do not find a dimension bound at any value of d, so that one can only see these 
bounds when a large number of derivatives are considered. 11 It would be very interesting to 
see if pushing the numerics further and incorporating even more derivatives could lead to 
bounds at larger values of d. 

We can also consider bounds on the OPE coefficients of operators without making any 
assumptions about the spectrum. In this case we simply require that a(J"A.o) > for all 
A > 2, which is the SUSY unitarity bound for scalar operators with vanishing U(1)r charge. 
In Figure 5 we show the resulting bounds on |Ao | for scalar operators appearing in this 
OPE as a function of their dimension, at various values of d. This is a supersymmetric 
generalization of the bounds considered in [35] in non-supersymmetric theories. Here we 
see that the bounds become very strong as A is increased, and appear to approach zero 
asymptotically. On the other hand, there are still finite bounds at Ao = 2, which tells us 
that even the coefficients appearing in front of flavor symmetry currents cannot be too large. 
We will explore this in more detail in the next subsection. 

4.2 Flavor Currents 

When is charged under a flavor symmetry, Ward identities guarantee that flavor currents 
appear in the OPE 0x0*, and thus contribute non-trivially to the conformal block expansion 
of (00*00*). In this subsection, we review the relevant Ward identities for both general CFTs 
and superconformal theories, and compute these conformal block contributions. In the next 
section, we present bounds on these quantities. Although we will eventually specialize to the 
case of a single operator 0, let us first consider a collection of operators 0j transforming in 
some general representation under a flavor group G. 

Suppose G has generators T 1 and conserved currents J Ia . The flavor charges in radial 
quantization are given by integrating the radial component of J Ia over a three-sphere sur- 
rounding the origin, Q 1 = —i J dVtx ■ J 1 . These charges then act non-trivially on our 0j as 
Q 7 0i(O) = —{T I )l 0j(O). 12 Comparing this action with our expression for Q 1 , we see that the 

n As discussed in Appendix B, Wk has dimension tih±ll^ so that k = 6 corresponds to 21 derivatives. It 
may be that not all of these derivatives are important for obtaining a dimension bound, and one possible 
numerical optimization might involve using a subspace of Wq other than Wk for k < 6. 

12 A word about i's and — l's. The minus sign in the action of Q 1 on <f>i ensures that commutators of 
Q 7 's act correctly. The — i in the definition of Q 1 comes from Wick- rotation to Euclidean signature. This 
is easiest to see in the usual time slicing, where / d 3 xJ° — > —i J d 3 xJ^ under J° — > —iJ%- OPE coefficients 
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max(A $t $) 

A 

^^^^^^^^^^^^^A = 2d 

d 

1 1.025 1.05 1.075 1.1 1.125 1.15 

Figure 4: An upper bound on the dimension of (the lowest-dimension scalar appearing in 
the $ x OPE), as a function of d = dim$. Here, we have taken = 6. The bound appears 
to approach 2 as d — > 1, as expected. On the other hand, we do not find a dimension bound for 
cZ > 1.16. It is possible that pushing the numerics beyond k = 6 could lead to bounds in this region. 



max | Xq | for Z = 




Figure 5: An upper bound on the OPE coefficient |Ao 1 of a scalar primary Oq appearing in $ x 
as a function of Ao = dimOo, f° r d = dim^ = 1.05, 1.1, 1.25, 1.5 (with k = 5). Here we no longer 
assume Oq is the lowest-dimension scalar in 3> x <J>t. 



3.75 
3.5 
3.25 
3 

2.75 
2.5 
2.25 
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J x OPE must take the form 

where ". . . " represents other operators, and we have used that vo^S* 3 ) = 2tt 2 . 
Suppose the 0j's and J /a 's are normalized so that 

(<MxiM(x 2 )> = % and (4.4) 

X 12 

Together, Eqs. (4.3) and (4.4) give the three-point function 

?' r 2 ~ 2d 

(MziWWW) = -TT-^-t^Z^ (4.6) 

where T~ = (T 1 ) 3 ^^. Combining this with Eq. (4.5), we find that the conformal block 
corresponding to an exchange of flavor currents in the <fi x <fi* channel is given by 

x^xlH^Mj) IrjjT^Tig^ (general CFTs), (4.7) 

where Tu is the inverse of r IJ . 

In super conformal theories, flavor currents J Ia (x) are the 9a a 9 terms in scalar supermul- 
tiplets J ! {z) of dimension 2. Comparing Eq. (4.7) to the superconformal block Eq. (3.43) 
with A = 2 and I = 0, we see that flavor supermultiplets contribute to a four-point function 
of anti-/chiral superconformal primaries as 

xljxHiMiMj) ~ T ijTfiTjj £?2,o (SCFTs). (4.8) 

Although the coefficients t ij are incalculable in general CFTs, in superconformal theories 
they have a simple expression in terms of the U(1)r generator [48, 49]: 

T IJ = -3Tr (i?T 7 T J ), (4.9) 

where the trace stands for the coefficient of the U(1)r anomaly induced by weakly gauging 
the flavor currents J Ia . For those SCFTs which emerge from a weakly coupled UV theory, 
this can often be calculated via 't Hooft anomaly matching. 

that don't involve e-tensors are the same in Euclidean and Lorentzian signature, so we are free to compute 
them in either signature. 
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4.2.1 Flavor Bounds 



Consider now a single scalar primary </> = 1; normalized so that g x j = \. We can bound the 
flavor current contribution tjjT^jT^j using the same procedure described in Section 4.1, with 
slightly modified constraints on the linear functional a. First, we demand that a(F^j) > 
(or «(J"a,«) > in the supersymmetric case) for all pairs (A,/) obeying the relevant 
unitarity bound. In general, this is A > 1 when / = 0, and A > / + 2 otherwise, while 
in a supersymmetric theory, it is simply A > 1 + 2. We also require a(F 31 ) = 1 in the 
non-supersymmetric case and a^J^.o) = 1 in the supersymmetric case, since these are the 
conformal blocks whose coefficients we wish to study. Note that we are no longer making 
implicit assumptions about the spectrum of operators appearing in x 0*, so the resulting 
bounds hold in any unitary CFT with a charged scalar primary. 

An upper bound on tjjT^T^j as a function of d = dim^ is shown in Figure 6 for a 
general CFT, and Figure 7 for a super conformal theory. Both bounds are strongest when d 
is near 1, and become weaker as d increases. The supersymmetric bound is most stringent, 
requiring tijT^T'^ < 1.6 when d ~ 1. 

Let us pause for a moment to appreciate the non-trivial nature of these bounds. If we 
for example consider a global U(l) symmetry with charges Qi in an asymptotically free 
superconformal theory, then using Eq. (4.9) we are placing an upper bound on the quantity 

Q 2 



-3£i(fli-l)Q 



2 ' 



(4.10) 



where the sum runs over chiral superfields in the UV description, and we are considering a 
gauge- invariant operator with charge Q. First note that this quantity does not depend on 
the overall normalization of the U(l) charges, which is unphysical. Our bound immediately 
tells us that one cannot have a global U(l) symmetry that acts only on fields that have R 
very close to 1. In addition, in principle one could imagine a cancellation between terms in 
the denominator, since some fields may have R smaller than 1 and some may have R greater 
than 1. Our bound also tells us that an arbitrary cancellation between terms is not possible. 

These bounds are also potentially interesting in light of the AdS / CFT correspondence [50- 
52]. In this case t ij is directly mapped to the size of the coupling constants for the 
corresponding bulk gauge fields. In AdSs this mapping is given by [53] 

r IJ = 8tt 2 L((T 2 ) /j , (4.11) 
where L is the AdS length scale and the gauge coupling (g~ 2 ) IJ appears in the action as 

'1 



SAdS = d Xyf^g 



-(9~ 2 y j FrFj, u + ... 



(4.12) 



Our bound tells us that there is a fundamental obstruction to making the gauge couplings 
arbitrarily large in the presence of charged scalar bulk excitations corresponding to operators 
with dimension close to 1. It would be very interesting to explore this connection further in 
a controlled setting, and to see if there is any kind of bulk reasoning that could give rise to 
this bound. 
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max(r/jT 7 T J ) for a charged scalar 

tuT 1 T j 




Figure 6: An upper bound on tjjT^T^- which is —3 times the flavor current coefficient in the OPE 
(f> x (p* of a complex scalar with its conjugate, as a function of d = dim (p. Here k = 5. 



max(r/jT 7 T J ) for a chiral primary 




Figure 7: An upper bound on the flavor current coefficient tuT^T^j appearing in the OPE 0x0* 
of a chiral primary with its conjugate in an N = 1 SCFT, as a function of d = dim0. Here k = 5. 
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4.3 The Stress Tensor 



The stress tensor T ab also makes a non-trivial contribution to conformal block expansions. 
Let us now review the relevant Ward identities, and compute the coefficient of the stress 
tensor conformal block in the supersymmetric and non-supersymmetric case. In the following 
subsection we will present bounds on these contributions. 

The dilatation operator in radial quantization is given by D = (—i) 2 J dQx a Xi,T ab , where 
the integral is over a three-sphere surrounding the origin. Requiring the action D<ft(0) = 
d(p(0) then determines the OPE 

T ab ( X )m ~ -y^(^ fc -i^ 2 )^(°) + --- ( 4 - 13 ) 

The stress tensor is conventionally normalized as 

(T ab (x)T cd (0)) = ^ /aC(3;) f d(x) , (4.14) 

7T 4 X 8 

where we implicitly symmetrize and subtract traces in each pair of indices a, b and c, d on 
the right-hand side. The coefficient c is the central charge, which appears for example in the 
trace anomaly (T£) = y^(Weyl) 2 — ^^(Euler) of the theory on a curved background. In 
our conventions, a free real scalar has c = ^ while a free Weyl fermion has c = ^. Just as 
with the flavor current normalization r /J , there is an explicit formula for c in superconformal 
theories, 

c = — (9TrR 3 -BTrR), (4.15) 

where R is the U(1)r generator, and the traces stand for anomaly coefficients. For a free 
chiral superfield (R = 2/3), the above equation yields c = ^. 

Now combining Eqs. (4.13) and (4.14), we obtain the stress tensor conformal block 
contribution 

x?^<#*#*) ~ ^4,2, (general CFTs) (4.16) 



where we've assumed that is normalized to have (<j)(x)(j)*(0)) = x 
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In superconformal theories, the stress tensor is the 9a a 8 component of the supercurrent 
J^ a (z) [54], a supermultiplet with dimension 3 and spin 1. Comparing with the superconfor- 
mal block Eq. (3.43), we see that the supercurrent contribution is 

xlixliiW*^*) ~ -^03,1 (SCFTs) (4.17) 

Note that the lowest component of J a is the U(1)r current J^, which has a conformal 
block contribution dictated by Eq. (4.7) with T R = Comparing with Eq. (4.17), we find 
t rr _ i|c ^ j s indeed correct (see for example [55]). 
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min(c) for general CFT 

c 

0.01 r 




Figure 8: A lower bound on the central charge c in a general CFT containing a real scalar primary </>, 
as a function of d = dim^. Here, we show the optimal bound for various values of k, corresponding 
to restrictions of the search to different finite-dimensional subspaces Wfc C V* (Eq. 2.15). As k 
increases, the bound gets stronger, though it becomes more computationally intensive to obtain. 
The bounds possibly converge to the free scalar value c = 1/120 as d — >■ 1. 



4.3.1 Central Charge Bounds 

We can now produce bounds on the central charge c using the same procedure as for flavor 
currents in Section 4.2.1, but with the equality constraints modified to a{Fi^) = 1 in the non- 
supersymmetric case and 0(^-3,1) = 1 in the supersymmetric case. Note that the coefficients 
in Eqs. (4.16) and (4.17) are inversely proportional to c, so that an upper bound on conformal 
block coefficients implies a lower bound on the central charge c > f c (d), as a function of 
d = dim <p. 

We plot this bound for the case of a real scalar in Figure 8. We have included curves 
for different values of k (indexing the size of our finite-dimensional subspaces Wk) to show 
how the bound gets stronger as we widen the search space S fl Wfc. In particular, as k 
increases, the series of bounds c > fc k \d) appears to approach c > Cf rcc scalar at d = 1. Recall 
that precisely at d = 1, the operator is free and decouples from the rest of the theory, 
contributing exactly Cf ree sca iar to c. We conjecture that / c (1) — > Cf rce sca i ar as k — > 00, namely 
that the optimal bound at d — 1 can be achieved with these methods. 

In Figure 9, we plot a lower bound on the central charge c > /f USY (<i), in any supercon- 
formal theory containing a scalar chiral primary of dimension d. Incorporating constraints 
from superconformal symmetry into the crossing relations certainly gives a stronger bound 
than in the case of a real scalar. However, we do not have (1) ~ Cf ree chirai superfieia, 
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Figure 9: A lower bound on the central charge c of a superconformal theory containing a scalar 
chiral primary <f>, as a function of d = dim <j>. Here we have taken k = 5. Note that the assumption of 
supersymmetry allows us to strengthen the bound significantly from the real scalar case (Figure 8). 
However, the bound does not appear to approach the free chiral superfield value c = 1/24 near 
d = 1. This is likely because we are only using the partial crossing relation Eq. (2.22) instead of 
the full information in Eq. (2.20). 

possibly reflecting the fact that we are using only the partial crossing relation Eq. (4.1). 

Unfortunately these central charge bounds are not "additive." That is, they are not 
stronger in the presence of multiple degrees of freedom, unless those degrees of freedom are 
completely decoupled from one another. Consider, for instance, a CFT with n real scalars 
transforming in the fundamental of an SO(n) flavor symmetry. If these scalars are decoupled, 
each with central charge c, then we can safely write 7f"^ ltheory = ^ i T- lb , where Tf h is the 
stress tensor in the z-th decoupled sector, and compute Cf u n theory = nc. In this case, a lower 
bound on c translates trivially into a much stronger lower bound on c^u theory However, 
now suppose our scalars are weakly interacting. We no longer have separate conserved stress 
tensors T" 6 , but rather a single stress tensor T ab which is a singlet under SO(n), along 
with a non-conserved spin-2 operator T^ kc transforming in a traceless symmetric tensor of 
SO(n). We can no longer say what the dimension of T^ e is, since it's no longer protected 
by a conservation law, and consequently we cannot straightforwardly include it in a linear 
program. 

While our bounds on c are perhaps somewhat weak, they are still highly non-trivial. 
In superconformal theories, for example, there is no a priori reason to think that different 
contributions to the central charge in Eq. (4.15) cannot cancel each other to a high degree. 
However, our bound says that this is not possible if the theory contains a scalar chiral 
primary operator of low dimension. 

In the context of AdS/CFT, c can be related to the bulk Planck scale Mp and AdS 
length scale L as c ~ ii 2 L 3 Mp [56]. Our bound then suggests that there is a fundamental 
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obstruction to making quantum gravity on A0IS5 arbitrarily strong in the presence of bulk 
scalar excitations corresponding to operators with dimension close to 1. It would be very 
interesting to make this more precise, and to understand the origin of the bound from the 
bulk perspective. 

Finally, we note that a somewhat similar bound on the central charge was derived in [36] in 
the context of 2D conformal field theories. There, an inequality relating c to the dimension of 
the lowest-dimension primary operator was derived through the use of unitarity and modular 
invariance, the latter of which is not available in 4 dimensions. 13 It would be interesting to 
derive bounds on c in 2-dimensional theories using the present techniques and see how the 
results compare to those of [36]. 

5 Comparison to Known Theories 

Since the SUSY flavor and central charge bounds derived in the previous section apply to 
quantities that are computable via 't Hooft anomaly matching using Eqs. (4.9) and (4.15), 
we can check whether they are satisfied in theories believed to flow to superconformal fixed 
points. Doing so requires knowing how the U(1)r subgroup of the superconformal group acts 
on our theory. In simple cases this action is uniquely determined by symmetry considerations, 
but in general it must be determined using a- maximization [57]. This requires knowing the 
full set of IR flavor symmetries that can mix with U(1)r, and in many cases accidental 
symmetries can arise that are not apparent from the UV description of the theory (for a nice 
discussion see [58]). In practice, one can sometimes identify the emergence of such accidental 
symmetries through apparent violations of the unitarity bound for the dimension of chiral 
primary operators, or by using a Seiberg dual description [59] if one is available. For a 
number of examples of such analyses, see e.g. [60-70]. In principle, apparent violation of our 
bounds could provide additional evidence for the emergence of such accidental symmetries, 
or even the absence of a superconformal fixed point. 

In this section, we will develop some intuition for the strength of our bounds by applying 
them to some simple superconformal theories, namely SU(iV c ) SQCD in the conformal 
window, and SQCD with an adjoint X and superpotential TrX 3 . In both cases, we find 
that the bounds are most interesting at small N c , although they are easily satisfied for all 
N c > 2 throughout the conformal window. While there are of course many other theories 
that can be checked, we will leave a more comprehensive survey of M = 1 theories to future 
work. 

5.1 SQCD 

Let us start by considering SU(N C ) SUSY QCD with Nf vector-like flavors Q l and Qj. For 
§iV c < Nf < 3N C the theory is believed to flow to an interacting conformal fixed point [59]. 

13 Or rather, modular invariance of the partition function on T = T 3 x S 1 doesn't lead to simple statements 
about the spectrum of operators, which is given by quantizing the theory on S 3 x time and not T 3 x time. 
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The anomaly-free global symmetries are SU(Nf)L x SU(Nf)u x U(1)b X Z7"(1) j r, with 
transforming as (Nf, 1,1,1 — jf) and Q? transforming as (l,Nf,— 1,1 — ^-). The ring of 
gauge-invariant chiral primary operators is generated by mesons M~ = Q l Qj and baryons 
B h...i Nc = Qh Q i Nc and ^ =Q~ i ... Q~ Nc . 

Mesons have dimension = 3(1 — N c /Nf) which can be close to 1 near the lower end 
of the conformal window Nf ~ \N C . Thus, they are good candidates for operators on which 
to check our bounds. Concretely, we may pick out a single component M\, and consider the 
constraints that crossing symmetry imposes on the superconformal block decomposition of 
the four-point function (MlM^M}M^). In particular, the Ml x M\ ] OPE contains both 
SU{Nf)i and SU(Nf)n flavor currents, along with the supercurrent J a . 

-U 



Let us focus on SU(A^j)^ and first compute r . We can work in the UV using 't Hooft 



anomaly matching. The fermions contained in Q have -R-charge — j^, so we find 



2 



T 



3 N 

13 - 3Tr(i?T 7 T J ) = —±8 IJ , (5.1) 



2N 



f 



where the generators T 1 of SU(iV/)i, have the usual normalization Tr(T 7 T J ) = \5 IJ in the 
fundamental representation, and the extra N c factor comes from summing over colors. Thus, 
we obtain 

rM = |i„T M = ^i, (5.2) 

where we have used the contraction 5ijT~T~ = ^{SijSjj— j^jdnSjj). Now, the meson Ml also 
gets an equal contribution from SU(Nf)n, so the total contribution from flavor currents is 
2(N f - l)/3iV c 2 . Note that this scales as ~ 1/N C for large N c and fixed N f /N c , so the bound 
is mainly interesting for small iV c theories. However, it is readily verified that all values of 
Nf and N c within the conformal window satisfy the bound given in Fig. 7. For example, 
taking N c = 2 and Nf = 4 we have du = 1-5 and a coefficient of .5, whereas the bound tells 
us that the coefficient cannot be larger than ~ 6. Similarly, taking N c = 3 and Nf = 5 we 
have du = 1-2 and a coefficient of .3, whereas the bound is ~ 2. Thus, while these theories 
are a factor of a few away, they do not come very close to saturating the bound. 

Finally, let us compare the central charge to the bound given in Fig. 9. We can calculate 

1 1 / TV 4 \ 

c = - (9Tri? 3 - 5Tri?) = - ( 7iV c 2 - - 2 j . (5.3) 

Note that this grows like ~ N^ for fixed Nf/N c , so the bound is again most likely to be 
interesting for small N c theories. However, it is also interesting that there are contributions 
to c with opposite signs, so in principle there could have been a cancellation. This occurs at 
Nf ~ 3N^ j \f7N^ — 2, which is always outside the conformal window. On the other hand, 
for all values of N c and Nf inside the conformal window c is greater than 1, and hence easily 
satisfies the bound in Fig. 9. 
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5.2 SQCD with an Adjoint 



Let us next consider SU(N C ) SUSY QCD with Nf flavors and an adjoint X. For simplicity, 
we focus on the theory with superpotential W = TrX 3 , which was studied in detail in [71]. 14 
This theory is believed to flow to an interacting fixed point for |iV c < Nf < 2N C . The 
anomaly-free global symmetries are SU(Nf) L x SU(Nf) R x U(1)b x U(1)r, with Q l trans- 
forming as (Nf, 1,1, 1 — f^-), Qi transforming as (1, Nf, — 1, 1 — |^), and X transforming 

as (1, 1, 0, |). Here we will focus on the chiral primary "meson" operators Mj = Q l Qj and 

N~ = Q l XQj. The theory also contains the chiral operator TrX 2 as well as baryons built 

out of products of Q % and XQ l , and anti-baryons built out of products of Qi and XQ^. 
Note that when |iV c < Nf < N c , Mj possesses a U(1)r charge that appears to violate the 
unitarity bound, and the interpretation in this case is that this operator has decoupled from 
the theory and become a free field. 

Let us begin with the case of Nf > N c . The dimension of Mj is given by du = 3 — 
and approaches 1 for Nf ~ N c . We will consider the bounds arising from crossing symmetry 
of the four-point function (M\M\ M\M^) . Each flavor group SU(Nf)i,.R has 

r" = ft"- (5.4) 
Since M is a flavor bifundamental, the flavor current conformal block contribution is 

(ruT^T^+frjT^Ttin = ^Lzl. (5.5) 



c 




Again this scales as ~ 1/N C , and there is no violation of our bound for any choice of N c < 
Nf < 2N C . Further, when Nf > N c the central charge is given by 

-4J. (5.6) 

Again, for N c < Nf < 2N C we always have c > 1 and are unable to approach the bound. 

Next let us consider the range |A^ C < Nf < N c . In this case the meson M becomes 
a free field and decouples from the rest of the theory. In the dual magnetic description 
this is simply described by the superpotential coupling involving M flowing to zero rather 
than a fixed point value [71]. In the present description we may equivalently describe this 
situation by adding to the theory a superpotential Wlm = L{(M~ — Q l Qj), containing new 
gauge-singlet fields L and M [63]. When Nf > N c , L and M are massive and can simply be 
integrated out, with the L equation of motion setting M = QQ in the chiral ring. However, 
when Nf < N c the "mass term" LM flows to zero and M is no longer interacting. We are 
left with a single new interacting field L whose equation of motion now sets QQ to zero in 



14 Note that one can straightforwardly generalize this discussion to the case of W = TrX k+1 or a vanishing 
superpotential using the results of [60, 72, 73]. 
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the chiral ring, thus avoiding the unitarity constraint. One must then include both L and 
M when computing r and c via anomaly matching. Of course, we already know that the 
central charge is at least as large as the contribution from M, so we cannot learn anything 
new from this bound. Additionally, since M has decoupled from the theory we can take it 
to transform under separate flavor symmetries as compared to the interacting sector. 

Now we will investigate the flavor current constraints imposed by crossing symmetry of 
the four-point function (N^N^ NlN^) in this regime. For both flavor groups SU(Nf)L ; R 
(which no longer act on M, but do act on L), we have 

r IJ = (w f + l Nf - 2Nc ) SlJ > (5J) 

from which we obtain 

One can then verify that for all |iV c < Nf < N c the bound of Fig. 7 is satisfied. 

6 Conclusions 

Let us point out some possible directions for future research. First, the bounds obtained 
in this work can of course be improved with more refined numerical methods. In the case 
of operators transforming under global symmetries, it also seems possible that additional 
crossing constraints can be incorporated that were not utilized in the present study. It 
would be interesting to see if doing so could lead to even stronger bounds on the dimensions 
of non-chiral operators in superconformal theories, so that one could start probing more 
phenomenologically interesting scenarios such as those of [28, 29]. Another interesting appli- 
cation is to see if one can bound the lowest-dimension SU (2)-singlet operator in conformal 
technicolor models, as was extensively discussed in [33]. In an ideal world, by incorporating 
the full set of crossing constraints one could perhaps obtain bounds that scale with the size 
of global symmetry representations. For example, the central charge c should roughly reflect 
the number of degrees of freedom of a theory, so the presence of large flavor representations 
should signal larger c. 15 It would be nice to derive a bound that supports this intuition. 

Another goal is to try to find Af = 1 SCFTs that come closer to saturating the bounds 
on c and tjj^T 3 . If any violation of these bounds could be found, it would be evidence 
for the emergence of new accidental symmetries or perhaps the absence of a conformal fixed 
point altogether. It may also be possible to extend these results to Af = 2 theories where 
one could obtain even stronger bounds. 

Similar studies in different numbers of dimensions are also feasable. The extension to 
other even dimensions should be completely straightforward since the conformal blocks are 

15 Actually, the central charge a generally appears to be a better measure of the number of degrees of 
freedom [48, 49, 57]. However, c is constrained by c > |a in supersymmetric theories and c > in general 
[74], so c is large whenever a is large. 
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known, and it would for example be interesting to see what kind of central charge bounds 
can be obtained in 2D using the present methods. While closed-form expressions for the 
conformal blocks in odd dimensions are not currently available, it seems likely that one could 
still use recursion relations (as described in Appendix B) to efficiently evaluate conformal 
blocks and their derivatives. Bounds obtained in 3D might then be relevant for condensed 
matter systems. 

Finally, it would be fascinating to better understand the interpretation of these bounds 
in the context of the AdS/CFT correspondence. They suggest that there should be a 
fundamental limit to the strength of gravitational and gauge forces in the presence of light 
bulk excitations in AdSs. Since our bounds are most interesting for small N theories, it 
seems likely that one will have to go to a highly quantum regime in order to see these effects. 
Nevertheless, it would be interesting to see if there is any simple bulk reasoning that could 
shed light on the origin of these bounds. One might then hope that thinking about these 
issues could lead to a deeper understanding of the nature of quantum gravity. 
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A Conventions 

Our metric and spinor conventions are those of the rj a b = diag(— 1, +1, +1, +1) version of [75]. 
The Clifford relation is cr a afc + o^a a = —2rj ab , so that one can convert between vectors and 
bispinors as (x) a a — x a&aa an d ^ a = — |tr(a a x). These conventions agree with those of Wess 
and Bagger [76] and Osborn [45], with a single exception — the sign of a , which affects the 
coefficient of e abcd in products of cr's and it's. Specifically, we have 

a a a b a c = - V ab a c + r] ca a b - if c a a - ie abcd a d (this paper) (A.l) 
uV(T c = -r] ab a c + V ca a b - r] bc a a + ie abcd a d (W&B) (A.2) 

To convert between these conventions, one simply flips the sign of e abcd wherever it appears. 

For the M = 1 superconformal algebra SU(2, 2|1), we follow the conventions used in [77]; 
in particular we take bosonic generators to be anti-hermitian (that is, they differ from the 
usual definitions by a factor of i). This eliminates some factors of i from the commutation 
relations, somewhat simplifying the algebra in Section 3. 

Let us arrange the superconformal generators according to their dimensions and spins as 
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follows 

dim(X) 

+ 1 Pa 
+ 1/2 Qa Q a 

M a(3 D,R M. $ ^ 

1/2 S a S a 

-1 K a , 

where M al3 = (a ba e) a ^M ab and M a/3 = (a ba e) al3 M ab are self-dual and anti-self-dual rotation 
generators. 

The dilatation operator and U(1)r generator act as 

[D, X] = dim(X) A [R, X] = i r(X)X, (A.4) 

where X is any generator, dim(X) is given in the above table (A. 3), and r(X) is the i?-charge 
of X, given by +1 for X = S, Q, by — 1 for X — Q, S, and zero otherwise. The additional 
commutation relations of the conformal sub-algebra are given by 

[M ab , P c ] = P a r] bc - P b r] ac , [M ab , K c ] = K a r] bc - K b r] ac 
[M ab , M cd \ = r] bc M ad - r] ac M bd - r] bd M ac + r] ad M bc 

[K a ,P b ] = 2 Vab D-2M ab . (A.5) 

Rotation generators act on spinors as 

[M ab , X a ] = (a ab )JX^ [M^ X*} = (a ab ) & pC^, (A.6) 

where X a = S,Q and X = Q,S. Finally, the remaining non- vanishing commutation relations 
involving fermionic generators are 



{Qa,Q & } 


— —2ja a . P 


{S a , S a } — 




(A.7) 


[K ai Q a ] 


"FT/3 


[S a , P a ] = 




(A.8) 




= <%, 


[s a ,p a ] = 




(A.9) 




{S a , Q/3J = 2De a /3 


- 2M af} - 3iRe a p, 




(A.10) 




{T,Q P } = 2De^ 


- 2M^ + 3iRe^. 




(A.11) 



The relation between our conventions for the super-Poincare subalgebra, and those of 
Wess and Bagger is summarized by equating supergroup elements at each point (x, 0, 6) in 
superspace 



e x-P+9Q+eQ 



this paper 



e i(-x-P+9Q+eQ) 



(A.12) 

W&B 



In particular, component expansions of our superfields 0(x, 8, 6) = e x ' p+e Q +e ®O(0) are the 
same as component expansions in Wess and Bagger, with the only difference being an overall 
factor of % or — % in the action of super-Poincare generators. 
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B Implementation Details 



In this Appendix, we discuss some details of our implementation of linear programs for 
extracting bounds from crossing relations. We first manipulate the crossing relation into a 
useful form, and then discuss efficient methods for calculation. Finally, we summarize our 
choice of programs and parameters for generating the bounds in this paper. 



B.l Explicit Formulae for Linear Functionals 

Using the explicit expression (2.8), we can rewrite the crossing relations Eq. (2.9) and (2.22) 

as 

(z — z) [z — z) 

[(1 -z){\ -z)] d ~ {zz) d 



IA 



2 l 



kA+i(z)k A _i_ 2 (z) k A+t (l - z)k A -i_ 2 (l - z) 



{zzy- 1 [(l-z^l-z)}^ 1 

-(z<*z) 7 (b.i; 



where the left-hand side is the contribution of the unit operator, and the sum is over the 
appropriate spectrum of primaries appearing in the OPE. Note that in the charged-scalar 
case, the (—1)' factor in the conformal blocks cancels with the (—1)' in Eq. (2.22), so that 
odd-spin contributions are not qualitatively different from even-spin contributions. 

We could bring this into the form (4.2) by additionally dividing by the left-hand side and 
isolating the term in the sum corresponding to a particular operator Oq. We would then 
consider the space of linear functionals a : f(z,z) i— > J2 m +n<2k a mnd™d^ f (1/2, 1/2) applied 
to both sides. Note however that we get the same space of functionals if we do not first divide 
by the unit operator, since derivatives of a product are linear combinations of derivatives of 
the two factors. Thus, to implement the algorithm described in Section 2.2, we can simply 
compute derivatives at z = z = 1/2 of Eq. (B.l) as written. 

Because of symmetry under (z,z) -H- (1 — z, 1 — z) and antisymmetry under z -H- z, 
derivatives d™d^ at (1/2, 1/2) will vanish unless m ^ n and m + n is even. Further, it 
suffices to take m < n by symmetry. Thus, our are precisely defined as the space of real 
linear functionals 



a 



f(z,z)^ a mn d?dnf(l/2,l/2), (B.2) 



m+n<2fc 

m+n£2Z,m<n 



which has dimension tittll, \y e w [\\ write the coefficients a mn collectively as a vector a. 

From Eq. (B.l), we see that the building blocks of these functionals are derivatives of 
z 1 ~ d kp(z) at z — 1/2. These have an analytic expression in terms of hypergeometric functions 
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which we can derive by matching power series, 

r(/?)r(/3/2 + m) 2 z ^-d+m 



n=0 



E 

m=0 

oo 

E 



r(p/2) 2 r(fi + m) ml 

T(P)T(P/2 + m) 2 r(m + g/2 - (rf - 1) + 1) 
i=o r(/3/2) 2 r(/3 + m)V{m + /3/2 - (d - 1) + 1 - n) 

xf-J 1 , 7 , j . (B.3) 

V2/ n\m\ K ' 



Performing the m-summation finally determines the coefficients 

C n M = 2 n + (^-i)-/3/2 ^ r(/3/2 + 2 - d )^ m/2 + 2 _ d ^ jg/2) ^ + 2 _ d _ w> ^ 1/2) , 



T(/3/2 + 2-rf-n) 



Now 9^9^|(i/2,i/2) applied to Eq. (B.l) can be written 



(B.4) 



2 2, + n +m -l (ci _ 1)(ra _ m)r(1 _ (i) 2 | Ao |2 

r(2-d-n)r(2-d-m) ~ ^ ^A +M ^A-z-2, d C A+M G A _ z _ 2j(i j , (B.5) 

where we have assumed that m + n is even. Hence the objective function in our 

linear program is given by a y v • a, where v is a vector of values of the left-hand side of 
Eq. (B.5) for different m and n (depending on our choice of W&). Meanwhile, each constraint 
tt(-^A,z) > becomes u A) ; ■ a > 0, where u A) ; is a vector of values of 

- \ rin rim rim rin 1 It) a\ 

L°A+i,d°A-/-2,d — u A+l,d u A-l-2,d\ 



2 l 



for different m, n. 



B.2 Optimizations 

Before running each linear program, we must compile a list of u Aii ^ for all (Aj, U) in our 
choice of discretization D. We found in practice that simply evaluating the expression (B.4) 
for C™ d in Mathematica introduced a performance bottleneck. One possible remedy is to 
precompute values of C^ d for different m, (3,d, and then perform table lookups to compile 
the list of constraints u Aii / r However, because there are three parameters m,(3,d to scan 
over, this would require a lot of memory and some careful bookkeeping. 

An alternative approach uses the fact that z 1 ~ d k^{z) = up td (z) satisfies a simple differen- 
tial equation which implies a recursion relation for its derivatives. Using the hypergeometric 
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differential equation for 2 i ? i(/3/2, (3/2, (3, z), and conjugating the resulting differential opera- 
tor by 2 ; /3 / 2+1 ~ d , we find 

z\\ + z(2d{l -z) + z-2)— + {d- l)(d(l - z) + z - 2) - A A u M {z), 

(B.7) 

where \p = \fi(f3 — 2). Now taking n — 2 derivatives with respect to z and evaluating at 
z — 1/2, we find the recursion relation 

C£ d = 2(5 - 2d - n)^ 1 + 4 (2A/? + n(n - 3) - d 2 + 4d - l) C™" 2 

+ 8(n-2)(n + d-4) 2 C™- 3 . (B.8) 

This can be iterated to give 

C n M = P n (f3,d)2 d -%(l/2) + Q n (f3,d)2 d -%(l/2), (B.9) 

where P n and Q n are polynomials in /3 and d. Now Eq. (B.9) can be made extremely 
computationally efficient. We first determine P n and Q n for all < n < 2k using (B.8). 
Additionally, we precompute a table of kp(l/2) and k'Jl/2) for different (3 values. This 
reduces the evaluation of C 1 ^ d to simple polynomial and exponential evaluation, along with 
two table lookups. 



B.3 Programs and Parameters 

Here, we give an account of the programs and parameters used to generate the bounds in 
Section 4. In each linear program, we take a discretization of the form 

D = {(A min + ne,Z) : n = 0, . . . , iV and / = 0, 2, . . . , L}, (B.10) 

where A m j n depends on the problem at hand, as discussed in the text. In addition to the 
parameters e, N, and L, one must also pick a subspace V\4 C V*. Our choices in this paper, 
along with the resolution of our plots are as follows: 
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4 


max(A <I) t$) 


50 


25 


0.02 


6 


Sd = 0.0025 


5 


max Xo for l = 


30 


25 


0.02 


5 


5A = 0.025 


6 


ma.x(rjjT I T J ) for a charged scalar 


20 


25 


0.05 


5 


8d = 0.05 


7 


ma.x(TijT I T J ) for a chiral primary 


20 


25 


0.05 


5 


6d = 0.05 


8 


min(c) for a general CFT 


80 


25 


0.02 


3,4,5,6 


6d = 0.01 


9 


min(c) for a SOFT 


20 


25 


0.05 


5 


5d = 0.01 



We have chosen iVe and L large enough so that the optimal linear combination satisfies 
®*(F&j) > asymptotically as A, I — > oo. At any finite e > 0, one can expect violations of 
the constraints a(F&j) > of order e 2 at isolated A, I. Decreasing e reduces these effects, 
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but has a computational cost since the linear programming algorithm we use (the simplex 
algorithm) runs in 0(l/e 3 ) time (cubic in the number of constraints). In each case above, 
we have verified that changing e slightly does not appreciably affect the results, so that we 
believe our plots accurately reflect the e — > limit. The curves themselves were generated 
by computing points with the resolution specified above (dropping a small number of points 
where the linear program was not well-behaved) and plotting an interpolating function. 

We generated the input data for each linear program with Mathematica. For actually 
solving linear programs, we used the GNU Linear Programming Kit (glpk) [78], which 
seemed generally faster and less unpredictable than Mathematica's LinearProgramming 
routine. Most of our computations were run on Harvard's Odyssey cluster supported by 
the FAS Sciences Division Research Computing Group. 
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